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ABSTRACT 


The  objective  of  this  thesis  is  to  investigate  the  problem  of  the 
safe  service  life  prediction  of  graphite  composite  structures  such  as 
pressure  vessels  used  in  the  pilot  ejection  seats,  rocket  motor  cases  and 
space  shuttle  energy  storage  compartments. 

The  basic  data  required  for  life  prediction  is  the  stress  rupture  life, 
i.e.,  the  composite  life  under  constant  load  and  it  has  to  be  recognized  that 
extensive  testing  over  long  periods  of  time  is  required  to  produce 
statistically  meaningful  data.  The  contribution  of  this  investigation  is 
focused  on  identifying  the  characterization  methodology  for  efficient  data 
generation  (to  minimize  cost  and  time),  the  appropriate  theoretical  models 
for  correlation  of  the  data  (to  translate  data  to  applications)  and  in 
evaluating  the  limitations  of  the  available  testing  (lower  bound  response 
speed)  and  data  acquisition  equipment  (upper  bound  recording  rate)  thereby 
enhancing  future  experimental  planning  and  testing. 
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INTRODUCTION 


A.         BACKGROUND 

Although  composites  have  been  employed  ever  since  materials  were 
first  used,  such  as  the  naturally  occuring  wood  or  the  man-made  mud  bricks 
reinforced  by  straws  and  have  been  widely  used  in  the  recent  50  years  in  the 
building  industry  (reinforced  concrete),  the  high  technology  of  the 
composites  has  evolved  in  the  aerospace  industry  only  in  the  last  25  years. 
The  first  strength  critical  application  for  modern  composites  were 
filament-wound  pressure  vessels  using  glass  fibers. 

What  has  highly  motivated  the  effort  of  composites  development 
has  been  the  goal  to  come  up  with  a  combination  of  properties  not 
achievable  by  any  of  the  constituent  materials  acting  alone.  Thus  a  solid 
could  be  fabricated  from  elemental  materials  which  by  themselves  could 
not  satisfy  particular  design  requirements  such  as  strength  and/or 
stiffness  to  weight  ratios  the  use  of  which  is  of  cruisial  importance  in  the 
aerospace  industry.  Another  highly  motivating  aspect  was  the  very  high 
reliability  of  the  composites  due  to  the  micro-redundancy  which. comes  as  a 
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result  of  the  fiber-matrix  load  sharing  that  in  fact  diffuses  the  local  fiber 
failure  sites.  And  it  was  this  materials  redundancy  that  finally  proved  thai 
composites  were  a  superior  application  for  monocoque  structures  such  as 
pressure  vessels  and  rocket  motor  cases. 

According  to  Tsai  [  Ref.  1  ]  the  results  of  this  effort  in  the  field  of 
the  Aerospace  Industry  came  sequentially  in  the  form  of  various  fuselage 
components  as  well  as  some  primary  control  surfaces  of  several  first  class 
fighters  such  as: 

1.  The  boron/epoxy  fuselage  section  and  horizontal  tail  of  the 
General  Dynamics  F-  I  1  1. 

2.  The  graphite/epoxy  fuselage  components  for  the  Northrop  F-5 
made  also  by  General  Dynamics. 

3.  The  limited  use  of  the  boron/epoxy  material  system  in  the  rudder 
of  the  Mc  Donnell-Douglas  F-4  which  has  accounted  for  a  35% 
weight  reduction  compared  to  its  Allurninum  counterpart. 

4.  The  boron/epoxy  horizontal  stabilizer  of  the  Grumann  F- 14. 

5.  The  composite  stabilator  of  the  Mc  Donnell-Douglas  F- 15. 

6.  The  graphite/epoxy  horizontal  and  vertical  stabilizers  of  the 
General  Dynamics  F- 16. 

On  the  other  hand  in  the  commercial  aerospace  industry  Boeing  was 
the  first  manufacturer  to  utilize  about  2  tons  of  composite  materials  in  the 
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floor  beams  and  all  the  control  surfaces  of  the  767  airliner  while  the  Beach 
Aircraft's  Starship  1  all-composite  aircraft  is  scheduled  to  start  its 
flight  test  program  before  the  end  of  1986. 

Apart  from  the  U.S  manufacturers  extensive  use  of  composites  has 
also  been  used  in  the  Dassault-Breguet's  Rafale,  the  Israel  Aircraft 
Industries'  Lavi  and  the  all-composite  fin  box  of  the  Airbus  Industrie 
A3 1 0-300,  an  impressive  structure  in  its  simplicity,  also  described  by  Tsai 
[  Ref.  1  I 

In  parallel  to  the  manufacturing  milestones  that  had  to  be  resolved 
once  the  goals  for  composite  compoments  were  set,  the  reliability  problems 
had  to  be  resolved  as  well.  Among  them  was  the  problem  of  safe  life 
prediction  of  the  composite  components  and  its  relation  with  the  composite 
strength.  The  significance  of  solving  this  problem  for  the  aerospace 
industry  becomes  more  obvious  when  one  focuses  on  the  construction  of 
space  and/or  missile  applications  where  repetition  for  experimental 
purposes  is  either  impossible  or  limited  to  numbers  often  less  than  3  due 
to  obvious  economic  constraints.  Consequently  one  has  to  realize  that  other 
methods  have  to  be  established,  methods  such  as  statistical  inference  and 
other  mathematical  techniques  which  will  enable  the  engineers  to  estimate 
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accurately  the  service  life  of  the  composite  component  or  system  under 
evaluation,  based  on  a  limited  number  of  experiments  which  is  acceptable 
from  the  cost  and  time  standpoints. 

From  physical  experiences  and  mathematical  modeling,  it  is  known 
that  the  average  macroscopic  strength  of  composites  depends  strongly  on 
the  number  of  micro-fiber  failure  sites  and  the  strength  of  specific 
composites  depends  on  the  clustering  of  the  fiber  failure  sites.  As  a  result, 
information  on  the  statistical  characteristics  of  fiber  filament  strength  is 
the  base  for  estimating  the  reliability  of  the  macro-strength  composite 
structures.  Specifically  the  statistical  lower  tail  distribution  of  the 
filament  strength  and  filament  life  are  required.  The  characterization  of  the 
lower  tail  behaviour  requires  testing  of  large  number  (thousands)  of 
samples.  In  the  case  of  strength  characterization,  large  number  of  samples 
is. a  matter  of  economics  whereas  in  the  case  life  characterization,  this  is  a 
matter  of  economics  and  time.  Aside  the  cost  constraints  which  are  always 
obvious,  the  reason  that  time  constraints  are  more  than  equally  serious  in 
this  case,  is  that  experimental  testing  of  fiber  filaments  could  last 
indefinitely  without  providing  adequate  data.  This  inherent  uncertainty  in 
life    characterization    can    only    be    overcome    by    accelerating    testing 
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methodologies.  In  this  investigation,  fiber  bundle  testing  is  chosen  to 
increase  the  number  of  samples.  However,  fiber  filament  testing  is  still 
required  to  serve  as  a  validating  tool  on  the  fiber  bundle  testing  results. 
Higher  stress  levels  are  used  to  accelerate  the  life  testing  in  order  to 
assure  that  failure  data  will  be  observed  within  the  existing  time 
constraints.  Under  stress-life  acceleration  circumstances  the  service  life 
prediction  could  be  made  possible  at  a  given  confidence  level  utilizing  the 
appropriate  statistical  methods,  given  that  the  fiber  bundles  under  testing 
had  already  been  able  to  withstand  either  a  prescribed  load  level  without 
any  failures,  or  a  higher  load  level  experiencing  a  known  number  of  filament 
failures  per  bundle.  The  results  of  this  investigation  are  relevant  to 
man-safe  applications  as  in  composite  pressure  vessels  used  in  the  pilot's 
ejection  seats,  the  liquid  oxygen  tank  in  the  jet  aircrafts,  the  space  shuttle 
tanks  for  energy  storage,  or  the  composite  rocket  motor  cradel  of  the 
guidance  system  of  various  missiles,  where  mostly  graphite  fibers  are 
used,  as  well  as  the  broad  class  of  composite  applications  where  the 
structure  is  subjected  to  sustained  tensile  loading. 

The  problem  of  accurately  determining  the  life  of  one  of  current  high 
strength  graphite  fiber  (Hercules  AS4)  has  already  been  started   at  the  MP. 5 
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by  Lt.Fred  D.Carozzo  and  the  results  of  his  thesis  [  Ref.  2  ]  identified  three 
major  milestones  that  have  to  be  eliminated  in  order  to  make  this 
characterization  methodology  operational: 


1.  To  avoid  the  stalling  of  the  tensile  testing  equipment  during  the 
loading  sequence,  because  stalling  appears  almost  exactly  as  an 
increase  in  broken  fibers  in  a  bundle  and  hence  can  produce  false 
interpetation. 

2.  To  explain  why  the  creep  data  recorded  during  the  creep-rupture 
phase  of  the  experiment  do  not  agree  with  the  expected  physics  of 
the  problem. 

3.  To  increase  the  recording  rate  of  the  data  acquisition  system 
used,  because  it  is  a  known  fact  that  much  more  fibers  per  bundle 
must  have  failed  during  the  loading  sequence  than  what  has 
already  been  recorded.  The  number  of  fibers  per  bundle  within  the 
bundle  which  failed  during  the  loading  seauence  is  needed  to 
calulate  the  conditional  probablity  for  the  AS4  fiber  bundle 
service  life. 


B.         GOALS 

The  objective  of  this  research  which  is  associated  with  the 
reliability  characterization  methodology  is  to  obtain  strength-life  data  for 
graphite  fibers.  The  final  target  is  to  present  these  data  in  such  a  way  that: 
"  Given  a  stress  level  one  could  predict  the  probability  of  safe  service  life 
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for  a  known  confidence  level  and  vise  versa  (that  is,  given  the  anticipated 
service  life  of  a  composite  application  specify  the  stress  level  that  should 
not  be  exceeded  during  sen/ice  in  order  to  achieve  the  given  safe-life  limit 
within  a  prescribed  confidence  level)  ".  With  this  ultimate  goal  in  mind, 
what  this  thesis  is  desired  to  contribute  to  is  four  fold: 


1.  To  present  a  general  lay  out  of  the  methodology  that  is  to  be  used 
as  a  guideline  to  reach  the  ultimate  goal  of  obtaining  strength-life 
data. 

2.  To  evaluate  what  is  the  minimum  loading  speed  that  the  existing 
testing  equipment  can  be  operated  with,  so  that  no  stall  will  occur 
during  the  loading  sequence  of  the  fiber  bundles  in  order  to  avoid 
interruptions  and  to  identify  the  load  dependency  of  the  minimum 
speed. 

3.  To  give  a  reasonable  explanation  on  why  the  creep  data  recorded 
during  the  creep-rupture  phase  of  the  experiment  do  not  agree 
with  the  expected  physics  of  the  problem. 

4.  To  study  the  limitations  of  the  existing  data  acquisition  system 
and  propose  solutions  to  increase  the  recording  rate  of  data  points 
during  the  short  time  loading  sequence. 


Consequently  the  projected  end-use  for  the  information  generated  bv 
this  investigation  can  be  specifically  tailored  to  the  applications  of  the 
Hercules    AS4   graphite    fiber   bundles    and    hence   the    upper    load    level 
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boundaries  for  specific  applications  can  therefore  be  established  such  that 
the  predetermined  service  life  of  the  composite  component  under 
evaluation  could  be  safely  achieved. 


2 


A 


APPROACH  TO  THE  PROBLEM 


A.    GENERAL 

It  is  intuitively  obvious  that  due  to  the  minute  transverse  dimensions 
of  the  composites  fibers,  huge  numbers  of  fiber  bundles  have  to  be  used  in 
the  construction  of  a  composite  structure,  regardless  of  its  absolute 
dimensions.  It  is  also  evident  that  the  type  of  constituents  to  be  used  in  a 
composite  structure  ought  to  have  been  experimentally  tested  in  advance 
and  demonstrated  that  they  comply  with  the  desired  standards  and/or 
specifications  concerning  their  strength  and-  service  life  prior  to  the 
beginning  of  the  manufacturing  process.  This  means  that  there  has  to  exist  a 
pilot  model  design  fabricated  under  known  parameters  prior  to  the 
commencment  of  the  manufacturing  process.  Furthermore,  since  the 
constituent  materials  used  for  the  model  evidently  differ  in  essence  from 
those  the  actual  manufacturing  will  start  from  (even  though  they  have  been 
fabricated  from  the  same  material  compositions  and  according  to  the  same 
procedures)  there  is  a  need  to  assure  the  reliability  of  the  production  line 
components.  The  answer  to  production  reliability  can  only  be  given  in  terms 
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of  confidence  levels  and  under  agreed  upon  conditions  between  the  Procuring 
Activity  and  the  Manufacturer.  It  is  also  likely  to  assume  that  within  the 
production  the  fiber  bundles  will  not  have  the  same  strength  or  the  same 
life.  Heuristicallly,  it  is  assumed  that  '*  high  strength  is  associated 
with  iong  life  whereas  low  strength  with  short 
fife  respectively  ".  Under  the  above  circumstances,  provided  it  is 
possible  to  identify  the  strongest  and  the  weakest  fiber  bundles  by  quality 
control  methods,  Figure  1  could  serve  as  a  good  starting  point  for 
understanding  how  strength  and  life  of  the  same  fibers  could  be  related, 
under  the  proviso  that  these  relations  have  to  be  determined  as  accurately 
as  possible    later  on.  Rosen  [  Ref.  3]    and  Phoenix  &  V/u  (  Ref.  4  ]  show  that 

both  the  strength  (at  t  =  tQ  where  tQ  is  a  very  short  finite  period  of  time), 

and  the  life  of  the  fiber  bundles    at  various  levels  of  strength    follow  the 
Weibull   distribution  which   is  usually   described  by   means   of   a  shape 

parameter  (aJ  and  a  scale  parameter  (j3j) 
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Figure  I,  WeiDuIl  Strength  to  weibull  Life  correspondence 
of  strong  and  weax  fiber  bundles  composed  of 
10  fiber  filaments. 
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B.         BASIC  Li NE  OF  METHODOLOGY 

Without  any  loss  of  generality    one  can  now    follow  the     thought 
experiment   described  below   for  a  fiber     bundle  composed   of  ten  (10) 

filament    fibers  identified  from  now  on  as  fj  for  i3 1,2,3,...,  10.  One  can  also 

assign  the  intrinsic  strengths  of  these  filament  fibers  to  be  (with  the 
implicit  assumption  that  the  numerical  values  of  the  strength  random 
variable  follow  the  Weibull  distribution)  as  follows: 


1.  One  filament  (assume  f  | )  of  strength  Pq]  =    1 0  lbs. 

2.  One  filament  (assume  f2)  of  strength  PS2  =    12  lbs. 

3.  Two  filaments  (assume  fj  and  f^)  of  strength  Ps^  =14  lbs. 

4.  Three  filaments  (assume  fc,  fg  and  fy)  of  strength  Ps^  =  1 6  lbs. 

5.  Two  filaments  (assume  fg  and  fg)  of  strength  P-^  =  18  lbs. 

6.  One  filament  (assume  f  jq)  of  strength  Ps6  =  20  lbs. 


If  one  desires  to  estimate  the  strength  of  this  bundle,  one  can  load 
the  bundle  to  different  strength  levels  (  PI,  P2,  P3,  etc.)  and  observe  the 
survival  of  the  filaments.  There  will  be  strength  levels  below  which  even 
the  weakest  fiber  of   the  ten  in  the  bundle  can  safely  withstand.  For  a 
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continuous  increase  in  stress  level    there  exists  be  a  load  level  (Pw<>  at 

which  the  weakest  fiber  will  finally  fail  and,  under  the  same  thought 
process,  there  must  be  another  load  level  at  which  the  next  weaker  fiber  in 

the  bundle  will  fail  (PW2^  ^  so  on,  until  all  the  fibers  in  the  bundle  will 

fail  and  hence  one  will  be  able  to  observe  the  failure  of  the  entire  bundle. 
Having  already  postulated  that  the  fiber's  strength  in  the  bundle  follows 
the  Weibull  model  distribution,  one  can  plot  the  probability  of  failure  for 

this  bundle  at  t  =  tQj  as  presented  in  Figure  2.    At  this  point  it  has  to  be 

noted  that  the  finite  time  t0  is  only  of  the  order  of  a  few  minutes  and 

therefore  negligible  if  compared  with  the  anticipated  duration  of  the 
service  life  and  therefore,  with  the  experiment  to  follow  for  the  life 
testing. 

A  very  important  observation  to  be  noted  is  that  the  fibers  in  the 
bundle  are  automatically  ranked  with  respect  to  their  strength,  which 
means  that  the  weaker  fibers  in  the  bundle  fail  first  whereas  the  stronger 
fibers  fail  last,  or  at  a  higher  load  level.  The  significance  of  this 
observation  is  what  actually  explains  the  big  advantage  of  the  fiber  bundle 
testing  versus  the  single     filament  r'iber   testing  and  can  become  more 
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Figure  2.     Homologous  correspondence  of  Weibull  Strength  at    t  =  to 

to  Weibull  Life  distnoutton   for  a  fiber  bundle  comoosea  of 
10  fiber  filaments. 
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obvious    if   one  thinks  of  bundles  composed  of  3000  up  to  10000  filaments 
(as  it  is  usually  the  case).  The  expediency   is  that  the   lower    tail    will 
always  be  the  first  to  be  observed.  As  it  has  already  been  stated  earlier, 

if  the  applied  load  level  does  not  exceed  for  example,  the  (P,v , )  level,  then 


no  failures  will  occur  until  t0.  However    if  one  were  to  leave  the  bundle 


loaded  at  a  load  level  (Pr. i<Pwi),  some  failures  could  be  observed  after  a 

long  period  of  time.  And  it  is  reasonable  to  think  that  under  this  load   P3, 

after  a  long  period  of  time  the  weaker  fiber  will  fail  first,  the  next  weaker 
fibers  will  fail  later,  and  so  on,  until  all  the  fibers  in  the  bundle  will  have 
failed,  when  it  would  also  be  possible  for  someone  to  observe  the  failure  of 
the  entire  bundle.  Assuming  again  a  Weibull  distribution  model  for  the  life 
of  the  fibers  bundle  as  per  Phoenix  &  Wu  [  Ref.  A  ]  one  can  plot  the 
probability    of  failure  of  this  very  bundle  in  time  or  in  other  words  in  life, 

given  that  no  failures  have  occured  during  the  loading  sequence  until  t0  * 


*  How  this  plot  could  be  obtained  will  be  explained  later  in  the 
Detailed  Method  of  Solution. 
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Usually  due  to  very  long  period  of  Lime  which  may  have  Ho  elapse  until 
the  first  failures  could  be  observed  it  is  customary  to  start  thinking  in 
Log  (  t  )  coordinates  as  far  as  the  abscissa  is  concerned  and  this  is 
the  way  the  figures  will  be  presented  during  this  thought  experiment. 

But  the  situation  described  so  far  is  merely  the  simplest,  and 
therefore  the  thought  experiment  has  to  be  carried  through  in  a  more 
realistic  way.  Such  a  real  world  application  could  be  one  in  which   the  load 

level  (P>„o)  and  this  load  level  could  have  been  maintained  to  be  constant  in 

vv  Z 

time  corresponding  to  a  stress  level  say  (Q9). 

Under  these  circumstances  what  should  happen  could  be  described  as 
follows: 

1.  During  the  loading  sequence  until  tQ  the  fibers  f  |  (Pw?  >  Pw  , )  and 
fo  (P,«o  >  Pr-^  )  will  fail,  whereas  the  remaining  (f-?  -f <...-.)  will 

survive. 

2.  Dunna  the  time  that  the  constant  load    P.A,n    is  maintained,  the 

3  W  Z  • 

rest  of  the  failures  have  to  occur  in  such  a  way  that  the  filament 
fibers  f-j  and  f^will  fail  earlier,  to  be  followed  by  the  failures  of 

filament  fibers  \'<=,  r>  and  f7,  the  failures  of  fQ  and  f'Q  and  finally 

the  failure  of  the  last  one  fjQ  which  completes  the  failure  of 

entire  bundle  . 


3° 


^ 


In  this  simplest  case  where  no  failures  had  occured  during  the 
loading  sequence  (due  to  the  low  magnitude  of  the  applied  load)  one  was  able 
to  infer  about  the  life  of  the  fiber  bundle  in  terms  of  the  Wetbul!  model 
based  on  the  magnitude  of  the  load  level  the  fiber  bundle  was  able  to 
withstand  without  failure.  Direct  statistical  inference  was  possible 
because  no  failures  had  occured  during  the  loading  sequence.  However  in  the 
real  world  applications  there  are  always  finite  number  of  failures  during 
the  loading  sequence  and  therefore  one  has  to  be  able  to  determine  the 
number  of  these  failures  and  their  respective  failure  stresses  in  order  to 
be  able  to  use  the  same  technique  and  infer  about  life.  In  other  words  life 
inference  is  possible  only  on  the  basis  of  the  conditional  probability 
that  failures  did  not  occur  during  the  loading  sequence.  And  in  order  to 
evaluate  this  condition  one  has  to  know  the  exact  number  of  the  fiber 
filaments  that  have  failed  during  the  loading  sequence,  so  that  the 
magnitude  of  the  load  level  that  would  have  been  achieved  if  no  failures  had 
occured  can  be  used  as  before.  This  requires  an  accurate  data  acquisition 
system  that  has  to  be  used  in  order  to  be  capable  to  record  the  possible 
failures  during  the  loading  sequence. 
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Although   the   sequence   of    failures   described   so    far   is   what    is 
anticipated  in  accordance  with  intuition,  (i.e  whatever  strength  of  the  fiber 

filaments  was  not  consumed  in  testing  during  the  loading  seduence,  can  be 

used     for     longer     life    endurance    under     a    homologous     strength-life 

relationship),  this  might  not  be  the  case    and  therefore  one  has  to  find  out 

experimentally  whether -this  one  to  one  correspondence  exists  between 

failures  at    t  =  t0    and  failures  in  time  (  t  -  t,  for  i  =  a,b,c,d  )  for  every 

arbitrary  selected  load  level  or  in  other  words  whether  the  above 
mentioned  relationship  of  strength  and  life  is  homologous  or  stochastic 
(see  Figure  2). 

On  the  other  hand,  another  major  constraint  is  the  fact  that  for  all 
practical  purposes  it  is  not  possible  nor  necessary  to  wait  until  the  failure 
of  the  last  filament  fiber  because  this  might  occur  much  later  than  the 
structural  lifetime.  This  is  the  reason  why  the  engineers  only  have  to  wait 
until  the  first  few  failures  occur  and  then  terminate  the  experiments  on 
purpose.  The  background  of  such  a  descision  lies  on  the  fact  that  the 
engineering  interest  is  in  the  early  failures  of  the  structure  under 
evaluation  which  justifies  a  safe  sen/ice  life,  in  terms     of  the  thoucht 

experiment  described  so  far    if  one  assumes  that  the  PW2  load  level  would 
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be  achieved  during  the  loading  sequence  at  t  =  tQ,  and  hence  fiber  filaments 

f,   and  f-  would  have  failed  and  furthermore  the  load  level  Pwo  can  be 

managed  to  be  reached  in  time,  one  should  have  to  wait  at  least  until  the 

first  few  failures  (namely  those  of  fiber  filaments  f^  and  f4)  will  occur 

when  the  experiment  can  be  terminated.  If  the  same  thought  experiment  is 
earned  out  several  times  reaching  higher  load  levels  in  every  loading 
sequence  and  infered  life  is  plotted  every  time  as  a  Weibull  distribution 
based  on  the  load  level  that  would  have  been  achieved  if  no  failures  had 
happened  during  the  loading  sequence,  a  curve  resembling  the  familiar  5-N 
fatigue  curve  will  arise  as  shown  in  Figure  3. 

The  significance  of  coming  up  with  a  graph  such  as  that  of  Figure  3 
for  a  real  word  application  as  the  widely  used  graphite  fiber  AS-4  is  more 
than  obvious,  apart  from  the  fact  that  serves  as  a  basis  for  the  so  called 
Proof-Test  with  the  intent  of  determining  definitively  both  the  strength 
and  the  life  of  the  same  specimen  simoultaneously. 
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Figure  3.    Weibuil    Life  distributions  for  Tour  identical 
fiber  bundles  comoosed  of  10  fiber  filaments 
acnieved  as  a  result  of  different  load  levels 
reached  at  tne  end  of  the  loading  sequences. 
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C.         DETAILED  METHOD  OF  SOLUTION 

What  remained  to  be  clarified  is  the  method  of  representing  the 
probability  of  failure  of  the  fiber  bundle  in  time,  given  that  no  failures  have 

occured  during  the  loading  sequence  (0  <t  <tQ)  where  t0  is  a  short  period  of 

time  of  the  order  of  several  minutes. 

To  demonstrate  how  this  can  be  done  and  also  get  more  insight  in  the 
process  of  approach  to  the  overall  problem  we  will  continue  with  the 
thought  experiment  considering  four  (4  )  identical  fiber  bundles  like  the  one 
which  has  already  been  described.  If  one  starts  loading  them  up  to  four  (4) 

different  load  levels  namely   Ps]  =  10  lbs,   Ps2  =  12  lbs,   Ps-j  =  14  lbs  and 

PS4=  16  lbs  and  then  maintain  these  four  load  levels  in  time  the  following 
pictures  can  then  be  presented  for  each  one  of  these  fiber  bundles: 

1 .         CASE  A .  10  filament  bundle  under  PG  |  =  10  lbs  tensile  load 

Since  Psj  is   the  strength  of  the  weakest  fiber  filament  in  the 

bundle,  obviously  no  failures  will  occur  during  the  loading  sequence. 
However,  after  some  period  of  time  failures  will  start  occuring  as  time 
increases,  giving  the  picture  shown  in  Figure  4- la.  Observing  carefully  the 
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steps  of  this  curve  it  will  not  be  difficult  for  one  to  understand  that  each 
one  of  these  steps  stands  for  certal"  failures  of  fiber  filaments  in  the 
bundle.  And  the  way  this  stepped  curve  is  plotted  shows  that  the  fiber 
filaments  of  the  same  strength  fail  simoultaneously.  Smoothemng  the 
stepped  curve  obtained,  it  is  not  difficult  for  one  to  observe  that  the  new 
smoothened   curve   stands   for   the   Reliability  Cummulative   Distribution 

Function  (RCDF)  of  the  fiber  bundle  under  evaluation.  And  since   the  Failure 

Cummulative  Distribution  Function  (FCDF)  is  always  related  to  the  (RCDF) 
by  the  relationship: 


(111)  (FCDF)  =  l-(RCDF) 


it  is  evident  that  FCDF  for  this  bundle  can  be  obtained  by  plotting  the 
complement  image  of  the  RCDF  with  respect  to  the  vertical  (load  )  axis  as 
shown  in  Figure  4- lb.  It  has  to  be  noted  at  his  point  that  this  picture  is 
usually  presented  in  the  so  called  Weibull  coordinates  due  to  several 
reasons  such  as  the  following: 


a.         The  Weibull  model  is  used  to  simulate  both  the  strength  and  th 
life  of  the  fiber  bundles  ([  Ref.  3  ]  and  [  Ref.  4  ]). 


39 


When  Weibull  coordinates  are  used  for  life  the  abscisa  is 
presented  in  Log  (t)  which  is  consistent  with  what  has  been 
already  stated  on  the  same  aspect. 

In  Weibull  coordinates  the  FCDF  represents  itself  as  a  straight 
line  instead  of  a  sinusoidal  shaped  curve  as  when  physical 
coordinates  are  used. 


Therefore  it  is  customary  to  present  the  picture  of  Figure  4- lb 
in  the  Weibull  coordinates  form  as  shown  in  Figure  4- 1c.  Taking  now  into 
account  that  no  failures  have  occured  during  the  loading  sequence,  one  can 
now  plot  the  failures'  probability  density  function  (fpdf)  using  the 
relationship: 


(11.2)  (fpdf)  -  [  d(FCDF)/dtj 


This  failures'  probability  density  function  for  this  first  fiber  bundle  denoted 
by  fpdf  is  finally  being  presented  in  Figure  4-  Id. 

2.         CASE B:  10   filament  bundle  under  P^o  =12  lbs   tensile  load 

Since    Ps2  is  the  strength  level  that  is  higher  from  only  one 

fiber  filament  strength,  it  is  reasonable  one  to  anticipate  one  failure  lo 
occur  during  the  loading  sequence  and  if  this  load  level  is  maintained  in 
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time  some  more  failures  will  start  occuring  until  the  whole  fiber  bundle 
will  fail  giving  a  picture  as  the  one  shown  in  Figure  4-2a.  T'r:e  only 
difference  between  this  case  and  the  previous  one  is  that  the  stepped 
curve  obtained  this  time  stands  for  the  Reliability  Cummulative 
Distribution  Function  of  the  fiber  filaments  that  survived  the  load  level 

achieved  at  the  end  of  the  loading  sequence  (RCDF)S.  In  order  for  one  to  be 

able  to  use  the  same  technique  as  in  the  previous  case  to  come  up  with  the 
failures'  probability  density  function  it  is  first  necessary  to  evaluate  the 
load  level  that  would  have  been  achieved  during  the  loading  sequence  under 
the  condition  that  no  failures  had  uccured.  This  load  level  can  oe  easily 
avaluated  given  that  the  number  of  the  fiber  filaments  that  have  failed 
during  the  load  sequence  is  known,  if  one  then  denotes  this  load  level  by 

PSNF  it  is  clear  that: 


(11.3)  (  PSjNF/PSj  )  =  (  100/%  of  survived  filaments  ) 


fon  i  =  1,2,3,...,  10 


Hence  fon  the  case  under  evaluation  (  i  =  2  ),  where  P.o  =  12  lbs  and    the 

3£ 


A\ 


03 


■a 

o 


CD 


C 

o 


0J 

t_ 

zj 

Q. 

i 

CL 
O 

U 

a 


■a 

c 

ZJ 


*-<  en 

Z>  *-) 

5  E 

TO  — 


CD 


CD    q    JD 

=     "~     CM 


ZJ 

3= 


o 


o 


ZJ 

.a 


CM 

XT 

CD 

ZJ 

cm 


fsq/jQVOl 


percentage  of  the  survived  fiber  filaments  is  0.90  (  one  failure  out  of  ten  ) 

one  can  easily  find  that   PS2N^  =  '3.3  lbs.  Then  the  upper  left  portion  of 

the  existing  stepped  curve  can  be  extrapolated  to  the  left  to  account  for  the 

(RCDF)^p    that  is  the  Reliability  Cummulative  Distribution  Function  under 

the  condition  that  no  failures  had  occured  during  the  loading  sequence.  Using 
the  same  reasoning  as  before  one  can  now  obtain  the    FCDF  for  both  the 

survived  fiber  filaments  in  this  case  (FCDF)S  and  for  the  fiber  bunble  under 

the  condition  that  no  failures  had  occured  during  the   loading  sequence 

(FCDF)^p.  Of  course  what  is  actually  important  for  this    case  is  only  the 

latter  that  can  be  evaluated  by  equation  (11.1).  This  result  is  presented  in 
Figure  4-2b  in  physical  and  in  Figure  4-2c  in  Weibull  coordinates. 
Finally  one  can  plot  the  fpdf  for  this  case  using  equation  (11.2)  as  presented 
in  Figure  4-2d. 

3.         CASE  C    10   filament  bundle  under  Pg3=14  lbs  tensile  ioao 

This  case  is  exactly  similar  to  what  has  already  been  described 
with  the  exception  of  the  load  or  stress  level  that  would  nave  D v ■. - r; 
achieved  at  the  end  of  the  loading  sequence  if  no  failures  had    occured 
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(PM,rsj    =    17.5    lbs).   Therefore    the   whole   procedure    is   outlined   only 
schematically  as  shown  in  Figures  4-3a  ,  4-3b  ,  4-3c  and  4-3d 

4.  CASE D\  10  filament  bundle  under  P=4=  16  lbs  tensile  load 

Again  due  to  the  fact  that  no  conceptual  differences  exist  no 
matter  what  is  the  creep-rupture    load  level  the  procedure  is  outlined  only 

schematically  (  PNFS4  =  26.7  lbs  ),  as    shown  in  Figures  4-4a.  4-4b, 

4 -4c  and  4-4d 

5.  Results 

If  one  now  takes  the  results  of  Figures  4- Id,  4-2d,  4- 3d  and 
4-4d  and  plots  them  with  respect  to  the  corresponding  load  levels  achieved 
at  the  end  of  each  loading  sequence,  the  picture  shown  in  Figure  5  can  be 
obtained  in  which  the  Strength-Life     curves  for  the  extreme  lower  tail 

(5-L)]t  ,  the  median  (5-L)m,  and  the  extreme  upper  tail  (5-L)ut  of  the 

Wcibull  model  distribution  are  plotted.    A  close  observation  of  this  figure 
can  lead  to  the  following  worthy  of  comment  thoughts: 


Given  that  one  can  produce  accurately  several  points  on  the  5-L 
curves  the  exact  determination  of  the  life  distribution  is 
possible  at  any  desired  strength  (stress)  level  is  only  a  matter  of" 
the  availability  of  equipment  if  time  is  not  a  problem. 
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Figure   5.    Strength   -  Life   curves  for  the  extreme  lower  tail  (5-Ult 
median  (5-Um   ana  uooer  tail  (5-Uut,  based  on  the  results 
of  four  experiments  of  identical  fiber  bundles  composed  of 
10  fiber  filaments  under  different  Creep-Puoture   loads. 
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D.        The    higher    the  number  or  fiber  filaments  in  the  bundles    the 
hiqher  the  confidence  level  of  the  results  to  oe  obtained,  or  in 
otherwords  the  more  certain  the  Weibull  curves    will    be      and 
with    the   less  spread  (i.e.,   smaller  coefficient  of  variation). 

c.  The  Weibull  curves  are  obtained  on  the  basis  that  one  can  produce 
an  accurate  curve  FCDF  under  the  condition  that  no  failures  had 
occured  during  the  loading  sequence.  The  accuracy  of  the  plot 
mainly  depends  on  the  exact  evaluation  of  the  fiber  filaments  that 
had  failed  during  the  loading  sequence  and  on  the  assumption  that 
every  fiber  is  carrying  the  same  amount  of  lead,  assuming  that 
internal  friction  is  absent  known  as  the  "equal  load-sharing 
rule  ". 

d.  What  is  of  great  practical  importance  is  the  exact  determination 
of  the  extreme  lower  tail  of  the  Weibull  distributions  for  various 
load  (stress)  levels  because  this  is  the  key  factor  for  evaluating 
the  safe  service  life  of  any  application.  Another  very  serious  point 
in  discussing  the  significance  of  the  extreme  lower  tail  is  the 
fact  that  the  time  coordinate  is  given  in  Log  (t)  and  therefore 
a  small  error  of  0.01  order  of  magnitude  in  Log  (t)  is 
translated  in  more  than  0.1  order  of  magnitude  error  in  useful! 
service  life  which  in  any  ca-.o  is  typically  unacceptable,  it  is  nor 
difficult  to  observe  that  the  exact  determination  of  the  extreme 
lower  tail  is  based  on  how  accurately  one  can  determine  the  load 
that  could  have  been  achieved  if  no  failures  had  hapDened  during 
the  loading  sequence.  Therer'ore  one  can  understand  that  it  is  of 
cruisial  importance  to  obtain  exact  data  points  with  no  irrelevant 
indications  during  the  loading  sequence  and  in  order  to 
accomplish  this  task  two  aspects  have  to  be  thoroughly  examined: 

I.         How  to  avoid  the  irrelevant  indications  when  evaluating 
the  data  points  recorded  during  the  loading  sequence. 
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How  to  record  as  m^ny   failure  data  points  as  possible 
to  increase  the  precision  of  the  final  results. 


D.         STATISTICAL  MODEL  OF  FAILURE 

The  failure  of  a  bundle  composed  of  brittle  fibers  such  as  carbon  and 
glass  in  a  flexible  matrix  under  tensile  loads,  is  a  complex  statistical 
process  involving  the  scattered  failure  of  fibers  at  imperfection  sites  and 
hence  the  overloading  of  the  adjacent  fibers  at  these  sites,  as  well  as  the 
propagation  of  the  neighboring  fiber  failures  to  a  critical  size  according  to 
Rosen  [  Ref.  6  ] . 

The  mathematical  mode!  fon  the  failure  process  was  finst  presented 
by  Rosen  [  Ref.  6  ]  and  is  described  in  detail  by  Phoenix  [  Ref.  7  ]  .However 
for  purposes  of  completeness  and  in  onden  to  gain  some  insight  in  the 
features  of  the  bundle  failune  a  review  of  the  simplest  case  is  presented 
henein. 

Acconding  to  this  model  the  composite  maternal  is  considered  to  be 
composed  of  repeating  representative  elements,  the  bundles,  which  are 
actually  the  focus  of  this  work.  Each  bundle  can  conceptually  be  viewed  as  a 
close    planar    (or    three    dimensional)    arrangement    of    n    parallel    fiber 
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filaments  that  form  a  tube  as  shown  in  Figure  6,  which  is  loaded  in 
tension.  This  tube  structure  is  in  turn  conceptually  partitioned  into  a  series 
of  m  snort  sections  tnat  are  called  microDunaies.  Furtermore  the  lengtn  <5 
of  these  short  sections  is  called  the  ineffective  length  after  Rosen  [Ref 
6]  and  actually  represents  the  length  of  the  ovenload  region  of  the 
neighboring  fiber  filaments  to  a  fiber  failure  site.  In  other  words  the 
composite  consists  of  bundles  each  one  of  which  has  a  length  equal  to  m£ 

and  a  volume  of  nmS.  When  a  moderate  tensile  load  per  fiber  p  is  applied  to 
a  bundle,  fiber  filaments  fail  randomly  at  sites  where  the  strength  is  less 
than  p  and  typically  these  random  failures  are  spatially  quite  far  appart.  At 
the  location  of  these  failures  sites  the  corresponding  fiber  filaments    can 

not  carry  any  load  and  therefore  an  overload  of  magnitude  Krp  takes    place 

on  each  neighboring  fiber     filament.,   where     Kr>  1    is  defined   as   the 

load  concentation  factor  and  r  represents  the  number  of  consecutive 
failed  fiber  filaments  immediately  adjacent  to  a  sound  fiber  filament 
counting  on  both  sides.  For  this  very  simple  case  that  is  being  described 


here  Kf  is  assumed  to  be: 
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Figure    6a   Statistical   Tensile  Failure  Model    after  Posen.W .B.IRef.  51 
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Figure  6  b.      Bundle  of  Tiber  filaments  :n  a  matrix 
in  trie  form  of  a  planar  taoe.  The  oundle 
is  comooseo  of   m  microoundies  in  series 
eacn   with   n    Tiber  filaments.  Failure 
sites  are  localized  within  microoundies. 


(11.4)  Kr  =  1+r/2       for    r  =  1,2,3,. 


The  length  of  this  overload  region  which  has  already  been  defined  as  8  is 
typically  of  the  order  of  a  few  fiber  filament  diameters  and  both  Kr  and  5 
depend  on  the  fiber  and  the  matrix  material  properties  which  can  be 
estimated  by  elasticity  or  numerical  analysis.  Under  the  increased  load  !<rp 
additional  failures  will  appear  to  the  close  neighboring  fibers,  which  in 
turn  will  be  subjected  to  a  higher  overload  say    KRp  (Kpp>Krp>  1 ).  This 

process  continues  until  all  the  failure  occurences  form  a  front  of  a  critical 
size  k*.  If  this  critical  size  is  exceeded  a  catastrophic  crack  propagates 
throughout  the  bundle  leading  to  the  failure  of  the  entire  composite 
structure. 

In  order  to  tailor  this  model  to  the  thought  experiment  already 
described  two  cases  of  stress  history  of  the  composite  are  of  interest. 

The  first  one  is  a  constant  stress  history  expressed  as  : 


(11.5)  o-j(t)  =  Q  for   t^O 
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where  Q  is  the  constant  positive  stress  level.  This  constant  stress  history 
is  often  referred  to  in  the  literature  as  stress  -  rupture  or  creep  - 
rupture  .  Application  of  this  type  of  stress  history  is  schematically 
illustrated  at  Figure  4- la. 

The  second  one  is  a  linearly  increasing  stress  history  expressed  as: 


(11.6)  <r2(t)  =  Rt  for    tiO 


where  R  is  a  positive  rate  of  stress  increase.  This  type  of  stress  history  is 
assumed  to  account  for  the  early  failures  of  the  fiber  bundles  that  will  take 
place  during  the  loading  sequence  and  it  is  often  described  as  short  - 
term  strength.  Application  of  this  stress  history  can  be  viewed  by  looking 
at  Figures  4- 2a,  4-3a  and  4- 4a.  The  FCDF  for  the  failure  time  of  a  single 
fiber  filament  of  length  8  is  given  by  Phoenix  and  Wu  [  Ref.  4  ]  after 
Coleman  [Ref.  8]  by  the  following  relationship: 


(117)  F(t;or)  -  1-expK[QJtx(cr(s))ds]]    ,     for     tiO 
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where  <r(t]  ,  t  )  0     is  the  stress  history, 

k[x]  -  yxP  ,  x  )  0    is     the  power-Jaw  breakdown  rule 

p   as  well  as  ?  are  positive  constants, 

$(x]  -  orxP  ,  x  )  0  is  the  Weibull  shape  function  to 
compensate  for  the  commonly  observed  Weibull  behavior  of  the  fibers  and 
a   as  well  as  p  are  again  positive  constants. 

Experimental  evidence  [  Ref.  4  Suggests  that  1  0  <  p  <  80   and 
that  1/4  <  p  <  4  for  a  wide  range  of  fibers.  It  has  also  to  be  noted  that  a 


a  I 


is  proportional  to  the  fiber-element  length,  or   a  =  <3r0S     for  a  consta 


in 


a:0  corresponding  to  the  characteristic  fiber  length  8Q.  Considering  8  to  be 
the  characteristic  length,  since  it  is  natural  to  the  composite  long  single 
fiber  filaments  of  length  1  >>  8  will  have  or  taken  as  mar  where 
m  \/8  The  theoretical  justification  of  this  fiber  model  is  discussed  by 
Phoenix  [  Ref.  7  ]  in  terms  of  the  Kinetic  failure  of  idealized  molecular 
crystals  of  the  form  found  in  stiff  poivrnenc  filaments  such  as  Kevlar  More 
specifically    the  key  parameter    p     involves  an    approximation    to    the 
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potential    function    for    chain    scission    and    is    shown    to    follow    the 


re'3-.:onsniD: 


(118) 


p  =  VkT 


where 


T 


R 


s  the  absolute  temperature, 


s  Boltzrnann's  constant 


U0  is   expressed  in  units  of  activation   energy  but  in   the  case 


01 


chain  scission  is  only  about  40  %    of  the  enerqv  of  the  bond 


I  UULUI  - 


whereas  the  other  parameters  <x  ,  p  and   y  vary  with  temperature  in  very 

complex  ways. 

Based  on  the  model  and  the  stress  histories  of  interest   described  so 
far    the  behavior  of  single  fiber  filaments  and  fiber  composites  will  be 


presented  from  both  the  sthrength  and  life  standpoints  together  wit 


h  the 


basic  assumptions  that  are  inherent  to  every  case  and  finally  some 
summarizing  results  will  also  be  presented  to  serve  as  a  means  of  gaining 
the  ^ciessary  insight  for  the  experimentally  anticipated  results  according 


to  Phoenix  and  Wu  [  Ref.  4  ]  and  Phoenix  [  Ref.  5 


1.        Behavior  of  single  fiber  filaments 
3.         Short-term  Strength 

The  fiber  filament  strength  is  known  to  oe  closely 
approximated  by  the  Weibull  distribution  which  under  the  linearly 
increasing  stress  history  (Eq.(ll.6)j  can  be  written  in  terms  of  the  (FCDF)  - 
F*  as  follows: 


9)  F*(p)  =  l-exp{-p/p,(R)]P(P  +  1)}    ,  for   p*0 


with  a  shape  parameter  p(p  +  1)  and  a  scale  parameter  equal  to: 


(II  9a)  p^R)  =  m~,/lP(PH)1p5<R) 


where     p^R]  -  a-VWp  +  M  [(i  +  p)/r]1/(?  +  1)  r1/(p  +  1) 


s  the  scale 


parameter  for  unit  fiber  length  8  and  R   is  the  loading  rate.  Finally  the 
median  fiber  strength  can  be  shown  to  be  equal  to: 
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(II  9b)  p*i(R)    =  p1(R)fln(2)J,/^(P+1)] 


Unfortunately  the  accuracy  of  Eq.(ll.9)  in  modelling  tne 
true  FCDF  can  only  be  demonstrated  for  the  lower  tail  of  F*(p)  that  is  for 

p<<Pj  and  this  is  because  strengths  associated  with  the  middle  and  upper 

tail  are  not  typically  observed  at  laboratory  gauge  length  1  (of  the  order  of 
several  centimeters).  However  this  does  not  constitute  a  serious  problem 
since  the  results  show  that  the  knowledge  of  the  middle  and  upper  tails  are 
not  important  in  the  characterization  of  composite  strength.  As  it  has 
already  been  explained  what  is  really  of  practical  importance  is  the 
accurate  determination  of  the  extreme  lower  tail, 
b.         Stress-Rupture 

A  single  fiber  filament  of  length  ]  =  m3  under  the 
constant  stress  history  described  by  Eq.(ll.8)  has  a  lifetime  which  can  be 
described  by  the  following  formula  after  the  appropriate  reduction  of 
Eq.(H.7): 


(1110)  F(t)=  1-expHt/tj(Q)]P]     ,     for   t>0 
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with  shape  parameter  p  and  scale  parameter  given  by  the  formula: 


(11.10a)  t,(Q)  =  m   1;P  U(Q) 


]\,\dJ    -    III  '       LQ\ 


where     td(Q)   =   Q  P/(:f<*^P)   is   the  scale  parameter  for  unit   fiber 

filament  length  8    For  the  constant  stress  level  Q  the  median  time-to 
-failure  can  be  evaluated  as: 


(11.10b)  t*,(Q)    =  t1(Q)[ln(2)]l/P 


if  the  median  strength  is  known  at  the  stress  rate  R,  and  one  wishes  to 
calculate  the  median  time  to  failure  at  the  stress  level  Q,  the  equivalent 
median  lifetime     can  be  evaluated  by  substituting  Q  in  Eq.(ll.lOa)  with 

p*j(R)  of  Eq.  (1 1.9b)  equal  to: 


(II. 10c)  t^R.p)  =  p*,(R)/lR(p*1)] 
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2.        Behavior  of  composite  fiber  bundles 
a.         Short-term  Strength 

One  can  now   consider  a  single  microbundle  and  define 

as      Gn(p)     for  p^O  the  (FCDF)  of  its  strength  or  in  other  words  the 

probability  of  failure  of  a  stochastically  selected  microbundle  under  a 
nominal  load  per  fiber  filament  equal  to  p.    Similarly  one  can  also  define 

Hm  n(p)  for  p*0    the  FCDF  for  the  strength  of  a  bundle  consisting  of  fiber 

filaments  in  a  matrix.  Since  the  necessary  and  sufficient  condition  for  the 
survival  of  the  bundle  is  the  survival  of  each  one  of  the  rn  microbundles, 
the  probability  of  survival  for  the  whole  bundle  under  the  load    p    can  be 

written  as    [  1  -Gn(p)]m    Therefore  the  probability  of  failure  of  this  bundle 

under  the  load   np,  amounts  to  the  so  called  weakest  link   rule  and  can 

be  presented  as  follows: 


(II  ID  Hm  „<d)  =   1-[1-G„(D)lm 


Hm  n(p)  =  1-[I-Gn(p)]m     ,     for   p*0 


Another  basic  assumption  that  has  to  be  stated  here  is 
the  fact  that  each  one  of  the  surviving  fiber  filaments  carries  the  same 
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amount  of  load  a  statement  which  amounts  for  the  equal-load  sharing 

rule  after  Rosen  [  Ref.  6  ].  The  first  ones  to  obtain  accurate  results  for  this 
kind  of  a  model  were  Harlow  &  Phoenix  [  Refs.  9  &10  ]  but  only  for  small 
number  of  fiber  filaments  per  bundle  (i.e.,  n<14  ).  Later  the  same  authors 
[Refs.  11  &  12]  developed  a  powerful  recursion  analysis  for  treatment  of 
larger  bundles.  According  to  this  analysis  the  FCDF  of  the  bundle  strength 
can  be  quite  succesfully  approximated  by  using  the  weakest  link  rule  as 
follows: 


(1112)  Hmjl(p)  =  1-[1-W(p)]mn    ,    forpiO 


where      W(p)    is   a  very   complicated   function  depending   on    the    load 

concentration  factor  Kr  and  the  FCDF  of  the  fiber  filaments  of  length  8 

under  the  tensile  load  p.  The  first  one  to  develop  a  tractable  approximation 
for  W(p)  was  Smith  [  Ref.  13  ]  who  showed  by  means  of  asymptotic 
analysis  for  large  n  (i.e  mn  =  106)  that  Eq.  (11.12)  can  take  the 
following  form  if  one  also  replaces  the  constant  stress  history  Q  (described 
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in  this  case  by  the  constant  load  p)  by  the  linearly  increasing  stress 
history  Rt : 


(1113)      *Hmn(p)  ^  1-expHp/pc(R)]k*P(P+l)}  ,  for  p^O 


with  shape  parameter     k*p(p+1)     and     scale  parameter  given  by  the 
formula: 


(11.13a)  pc(R)  =  (mnr,/tk*P(P+,)]  P<5(R)(dk*),/(P+  *  > 


where  k*  represents  the  critical  failure  sequence  size  and  mathematically 
is  the  integer  solving  the  equation  : 


(II. 13b)  l/r(k+l)  <  pp/ln(mn)  <  1/  r(k) 


where  the  function    r(k)   is  defined  by  the  following  equation 


(11.13c)     r(k)  =  ln(  k_1Kk  )  -  ln(  j=,  ]"]  k  KM )    and 
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dk-r(pHr,/Pr(kpH)1/<kP>2(,-k)/(kP)Jai|nk<kj-n"P/k 


for    k  =  1,2,3,...  CM.  13d) 


As  before  the  median  strength  can  be  given  by  the  equation: 

(11.14)  P*C(R>  =  Pc(R)nn(2)]1/[R*P(P+1)1 

b.         Composite  Stress-Rupture 

Similarly  the  FCDF  for  the  lifetime  of  the  composite 
material  under  the  constant  stress  history  described  by  Eq.01.5)  has  a 
Weibull  approximation  of  the  form: 

(11.15)  *Hm>n(t)  ~  1-expHt/tc(Q)]k*P  ,    for  t*0 

with  shape  parameter  k*p  and  scale  parameter  given  by  the  equation: 


(II.  16)  tc(Q)  =  (mnr1/(k*P)(dk^)t5(Q) 


where  d^*  has  already  been  defined  by  Eq.(ll.13d)  for  k*  =  k  and    t^(Q)  is 
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defined  in  the  same  way  as  for  the  Eq.(II.IOa).  The  median  lifetime  of  the 
model  composite  can  therefore  be  given  by  the  following  expression: 


(11.17)  t*c(Q)  =  tc(Q)[ln(2)]1/(k*P) 


Using  the  reasoning  given  in  the  analysis  of  the  fiber 

n 
filament  behaviour  if  the  median  strength  p  C(R)  is  known  for  the  stress 

rate  R  the  median  lifetime    under  the  stress  level  Q  -  p  C(R)  can  be 
evaluated  as: 


(II  18)  t*  (R.p)  -  p*  (R)/[R(p*1)] 


It  has  to  be  noted  that  the  Weibull  approximations  (11.13)  and  (11.15) 
are  only  valid  when  0p  is  large ,   say  greater  than  6,  as  well  as,  that  the 

best  Weibull  fit  to  Hm  n(t)  may  produce  a  value  of  k*  which  is  not  an 

integer  usually  when  the  solution  for  k  yields  almost  an  equality  on  either 
side  of  Eq.(ll.13b). 

Summarizing  the  important  characteristics  of  the  model  described 
one  has  to  note  that: 

I.        Comparisons  of  the  median  lifetime  and  strength  of  a  single  fiber 
filament  and  a  fiber  composite  are  actually  meaningless. 
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2.  The  variability  in  strength  and  lifetime  for  a  composite  is  much 
less  than  that  for  a  single  fiber  filament  primarily  due  to  the 
critical  crack  size. 

3.  The  size  effect  is  a  lot  milder  for  the  composite  relative  to  that 
of  the  fiber  for  the  strength  and  drastically  lower  for  the  lifetime 
again  primarily  due  to  the  critical  crack  size. 
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THE  CONTRIBUTION  TO  THE  SOLUTION 


In  accordance  with  the  approach  stated  in  the  Introduction  the 
required  methodology  appropiate  to  be  used  as  a  quide  line  in  order  to 
achieve  the  ultimate  objective  of  obtaining  strength-life  data  been 
presented  in  the  of  the  problem  approach  section. 

For  the  equipment  used  by  Carozzo  [Ref.  2  :  Fig.  3.4]  shown  in 
Figure  7  and  from  his  evaluation  of  its  performance  (summarized  in 
Figures  B.  I  and  B.2a  of  his  thesis  ([Ref.  2])  shown  in  Figures  8a  and  8b  one 
can  observe  that  although  the  desired  target  was  the  acquisition  of  data  in 
hundredths  of  seconds  (i.e  6000  data/min)  the  results  suggested  that  a  rate 
of  only  3-4  data/min  was  achievable.  In  an  effort  to  acquire  more  data 
during  the  loading  phase,  the  testing  traverse  speed  was  lowered  to  the 
minimum.  This  slowest  speed  setting  of  the  IN5TR0N  testing  equipment 
(at  13  grams/sec)  caused  the  drive  motor  of  the  INSTRON  to  stall  erratically 
[Ref.  2  :  pp.  29-30].  In  addition,  during  the  the  creep  rupture  phase  of  the 
experiment  the  time  dependence  of  creep  was  physically  observed  by  the 
upward  movement    of  the  moving  cross  -  head  of  the  INSTRON.  However 
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Figure    7      :    Block  Diagram  of  the  initial  Experimental  Set-Up 
after  Carozzo  [Ref.  21 
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Figure    8a    :    Load  -  Time  curves  for  Test  *1 
after  Carozzo  [Ref.  2]. 


67 


AS4  GRAPHITE,    SPOOL  019 


TIME    CMIN) 


Q 

i 

^ 

; 

; 

< 

"t  - 

: 

Q 
Q 

rt  . 

| 

i 

a. 

9— 

n 

- 

-  ,  ji  — 

■flU— 

■  ■■  &■ 

3— 

— .-%_ 

—a — 

n 

- 

— g 

3 . — 

— a_ 

.TV  — 

« 

■  — n — 

G r— 

— 5— 

3 

-> 

—3 

c 

; 

i: 

i 

U 

k 

15 

16 

17 

1{ 

! 

i< 

> 

2< 

5 

21 

TIME    (MIN) 


Figure    8b     :    Load  -  Time  curves  for  Test  *2 
after  Carozzo  [Ref.  2]. 
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the  data  recorded  to  support  this  observation  were  considered  completely 
unreliable  in  that  it  showed  that  the  cross-head  moved  up  and  down  which 
is  equivalent  to  an  entropy  decrease  ([Ref.  2  :  p.  31]).  These  observations 
were  what  actually  gave  rise  to  the  other  goals  of  this  work  which  are  to 
evaluate  the  limitations  of  the  testing  equipment  available  and  propose 
procedures  that  will  serve  in  eliminating  them  (i.e  explaining  why  the 
creep-rupture  data  were  unreliable,  avoiding  the  stall  ,  by-pass  the  problem 
of  inadequate  accuracy  of  the  extension  measurements  etc.),  as  well  as  to 
find  out  how  more  data  could  be  acquired  at  least  during  the  loading 
sequence  for  the  reasons  that  have  already  been  explained  in  the  detailed 
analysis  of  the  approach  to  the  overall  problem. 


A.         TESTING  EQUIPMENT  LIMITATIONS 

The  implementation  of  the  INSTRON  model  1000  Universal  Testing 
Equipment  requires  at  least  the  description  of  it's  operating  principles 
prior  to  investigating  it's  inherent  limitations. 

1.        Operating  Principles 

The  INSTRON    1000   is  an  120  ±  10V  AC  single  phase  50-60HZ 
motor  belt  driven  testing  device  capable  of     tensile  and  compressive 
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testing,  extension  measurements  and  the  testing  machine  has  been 
modified  to  perform  creep  monitoring  if  necessary.  The  load  measurements 
are  achieved  by  means  of  a  combined  Load  Transducer  (LT)/  Linear 
Variable  Differential  Transformer  (LVDT)  system,  installed  within  the 
moving  cross-head.  The  traverse  crosshead  measurement  is  achieved  by  a 
Shaft  Encoder  (SE)  and  the  creep  monitoring  by  means  of  a  Digital 
Processor/Controller  (DP/C)  that  can  compensate  for  the  cross-head 
adjustments  relative  to  a  predetermined  load  variation.  The  creep 
dispacement  was  specifically  recorded  by  an  additional  LVDT  that  was 
mounted  on  the  fixed  cross-head  of  the  IN5TRON  1000.  The  operating 
principle  of  the  LVDT  is  that  of  a  transformer  composed  of  a  housing  of 
the  transformer  coil  and  a  core.  The  insertion  of  a  ferrous  rod-core  alteres 
the  electro-magnetic  coupling  thereby  causing  a  change  in  the  output 
voltage.  Therefore  by  controlling  the  core  insertion  and  measuring  the 
output  voltage  one  can  evaluate  a  relationship  that  will  allow  for 
converting  output  voltage  to  displacement  and  vice  versa.  A  detailed 
description  of  the  major  units,  subsystems  and  enhancement  features  of 
the  testing  equipment  and  their  primary  functions  is  presented  in  Section  I 
of  Appendix  A . 
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2.  Determination  of  INSTRON  10QO  Limitations  by  Simulation 
The  strength  of  the  fiber  bundles  to  be  tested  is  of  the  50  lbs 
order  of  magnitude  and  the  maximum  utilization  of  INSTRON  1000  under 
the  existing  provisions  set-up  involves  testing  of  nine  (9)  bundles 
simoultaneously  (including  the  control  bundle).  Hence  the  load  level  of 
interest  is  of  the  450  lbs  order  of  magnitude,  an  estimate  that  dictates 
for  the  use  of  only  the  1000  lbs  LT  during  the  simulation.  Given  also  that 
the  graphite  fiber  bundles  ought  to  be  tested  in  tension,  because  the  area 
of  interest  is  ultimately  the  broad  class  of  composite  applications  where 
the  structure  is  subjected  to  sustained  tensile  loading,  as  well  as  that  the 
graphite  fiber  bundles  are  expensive  to  use  them  for  purposes  other  than 
the  actual  experiment  a  simulation  process  is  obviously  necessary.  During 
this  simulation  experiment,  the  above  mentioned  load  level  (450  lbs)  will 
have  to  be  easily  achieved  and/or  overcome  a  good  number  of  times  so  that 
several  conclusions  could  be  drawn  as  far  as  the  limitations  of  the 
equipment  are  concerned. 

The  objectives,  the  theoretical  model,  the  design  of  a  simple 
mechanical  set-up  that  had  to  be  used  with  the  INSTRON   1000  during  the 
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simulation   and  the  detailed  procedure  of  this  aside  task  are  presented  in 
Section  II  of  Appendix  A  . 

3.        Stability  of  the  LVDT 

The  background  for  checking  the  stability  of  the  LVDT  lies  on 
the  fact,  that  although  the  cross-head  of  the  INSTRON  must  have  moved 
during  the  creep-rupture  phase,  no  reasonable  output  had  been  recorded 
accordingly  and  the  output  (displacement)  was  measured  by  means  of  an 
LVDT  which  was  externally  mounted  on  the  INSTRON  1000. 

The  background,  the  objectives,  the  set  up  used  and  the 
detailed  procedure  of  this  separate  aside  task  are  presented  in  Section  III 
of  Appendix  A . 


B.         DATA  ACQUISITION  SYSTEM  LIMITATIONS 

One  of  the  major  problems  associated  with  the  graphite  bundle  life 
testing  as  presented  by  Carozzo  [Ref.  2  :  Fig.  Bl  &  B2a]  (Figures  8a  and  8b)  is 
the  slow  data  aquisition  rate  that  was  achieved  during  the  loading  sequence 
(0.04  readings/sec)  was  substantially  less  than  the  manufacturer's 
specification  of  about  2  readings/sec  [Ref.  15  :  p.  32].  For  completeness  and 
future  reference  purposes  the  experimental  set-up  used  by  Carozzo  during 
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the  initial  run  of  the  experiment  is  presented  in  Figure  9.  A  survey  of  the 
available  literature  on  this  subject  revealed  that  the  reading  rates  can  be 
influenced    by     many     factors. 

1.  Signal  environment  factors 

The  most  significant  factors  in  this  area  are  the  line  related 
factors  (the  smaller  the  loop  the  better  ,  60  Hz  line  frequency  gives  about 
17%  faster  reading  rates),  [Ref.  15  :  pp.  30,34]  &  [Ref.  16  :  p.  42],  the 
broadband  noise  (the  larger  the  loop,  the  more  the  noise  contrbutions  from 
the  various  accesories  in  the  loop),  [Ref.  15  :p.  30]  &  [Ref.  16  :  p.  42]  and  the 
thermal  gradients  (negligible  in  this  case  due  to  the  relatively  constant 
temperature  of  the  lab  environment),  [Ref.  15  :  p.  30]. 

2.  Desired  accuracy 

At  3.5  digits  of  resolution  (  N3)  the  ability  of  the  voltmeter  to 
acurately  measure  DC  voltages  in  the  presence  of  AC  voltages  at  power  line 
frequencies  is  expressed  by  the  Normal  Mode  Rejection  (NMR),  [Ref.  15  :  pp. 
30,  43],  which  in  this  case  is  0  db  NMR.  The  smaller  the  number  of  digits 
of  resolution  the  smaller  the  "Integration  Time"  which  obviously  determines 
the  reading  rate  and  hence  for  (N3)  the  reading  rate  is  faster.  Similarly  at 
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Figure  9  :  Block  diagram  of  the  experimental  set  -  up 
used  during  the  initial  run  of  the  experiment 
involving  one  IN5TR0N  1000  machine. 
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N4  (59  db  NMR)  an  intermediate  reading  rate  is  achieved  and  finally  at  N5 
(80  db  NMR)  the  reding  rate  is  slow 
3.        Type  of  measurement 

a.        Resistance  measurements  [Ref.  15  :pp.  30-32]. 


(1)  At  30  Mohm  (  200  msec  settling  time  )  the  recording  rate  is  1.5 
readings/ sec. 

(2)  At    3    Mohm  (  20  msec  settling  time  )  the  recording  rate  is  2.1 
readings/sec . 


b.        Voltage  measurements. 


( 1 )  The  DC  readings  are  5  to  50  times  faster  than  the  AC  ones  [Ref.  1 5 
:p.  32]. 

(2)  Positive  voltage  measurements  increase  the  reading  rate.  Possible 
problems  can  be  raised  by  the  LVDTs  whose  outputs  can  be 
negative  dependent  on  their  null  point  with  respect  to  the 
measurements  under  evaluation. 

(3)  AUTOZERO  is  a  function  that  allows  the  operator  to  enable  or 
disable  the  internal  zero  correction  of  the  voltmeter  of  the 
HP-3421A.  The  default  condition  is  ON  (  used  by  the  set-up  of 
Figure  9  )  but  turning  AUTOZERO  OFF  could  substantially 
increasthe  reading  rate  but  with  a  tradeoff  of  long  term  stability 
and  accuracy  [Ref.  15  :  p.  41].  However  no  problems  of  this  sort  are 
of  interest  during  the  relatively  short  period  of  the  loading 
sequence  which  is  of  the  order  of  a  few  minutes  (see  Figures  8a 
and  8b). 
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nherent  to  the  HP-3421A  features 


Calibration  has  to  be  performed  after  the  installation  of  the 
one  44462A  /  10  channel  Multiplexer  (option  020)  but  it  is  of  interest  only 
when  absolute  values  are  necessary,  however  this  is  not  this  case. 

AUTORANGING  is  a  function  meaning  that  signals  with  a  wide 
dynamic  range  such  as  switching  from  channel  to  channel  to  take  readings 
are  made  more  quickly  by  utilizing  it  [Ref.  15  :  p.  44].  But  this  function  is 
made  primarily  for  taking  readings  in  the  N4  mode  and  not  in  the  N5  mode 
which  has  to  be  used  due  to  the  HP-41CX  involvementfsee  Figure  9). 

When  operating  with  the  HEWLETT  PACKARD- INTERFACE  LOOP 
(HP-IL)  most  3421 A  commands  will  "  Hold  -  Up  "  the  computer  until  all 
readings  have  been  taken.  These  problems  are  known  as  information 
transfer  ones.  Recording  time  by  the  HP-41CX  when  the  order  of  magnitude 
of  the  3421 A  capacity  is  about  2  readings  /  sec  and  the  target  is  to  take  say 
100  readings  using  one  44462A  10  channel  multiplexer  assembly  50  seconds 
are  required  for  the  100  readings.  The  computer  will  not  be  able  to  perform 
any  operation  during  the  50  seconds  while  the  3421 A  is  taking  the 
measurements. 
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The  3421 A  commands  that  do  not  hold-up  the  computer  are  the 
so  called  advanced  commands  of  digital  monitoring  (MN,  MH  and  ML)  the 
totalize  function  (TOT)  and  the  digital  trigger  (DT)  commands.  However  if 
any  of  the  digital  monitor  modes  (MN  ,MH  ,ML  or  DT)  are  in  effect  in  the 
3421A  any  communication  through  the  HP-IL  will  cause  the  mode  to  be 
aborted.  In  other  words  if  the  3421 A  is  waiting  for  a  digital  trigger  and  the 
computer  or  another  instrument  in  the  loop  sends  any  commands  or  data 
through  the  interface,  the  3421 A  will  abort  the  digital  monitor  and  will  no 
longer  respond. 

On  the  other  hand,  according  to  Carozzo  [Ref.  2  :  p.25]  the  basic 
data  logger  software  program  provided  by  the  the  HP-44468A  Data 
Acquisition  Control  Pac  was  modified  to  permit  the  recording  of  time  in 
hundredths  of  seconds  and  edited  to  reduce  the  number  of  data  registers 
used.  However  recording  time  in  hundredths  of  seconds  is  50  times  faster 
the  maximum  capacity  of  the  HP- 3421 A  DACU  (2  readings/second).  In 
addition  to  this  the  system  timing  is  greatly  dependent  on  the  nature  of  the 
set-up  itself  and  especially  on  whether  the  printer  or  tape  cassette  (or 
both)  are  used  the  number  and  type  of  measurement  sequences,  User 
Functions  as  well  as  the  number  of  channels  in  each  sequence  [Ref.  16  .p.  42]. 
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Several  contacts  with  system  engineers  of  Hewlett-Packard  revealed  that 
the  reading  rates  of  the  HP-3421A  manual  for  the  case  of  the  experiment 
under  evaluation,  (Z 1 F 1  -2. 1 9/ 1 3. 1 8/22.22)  [Ref .  1 5  :  p.  32]  refer  only  to  the 
unit  itself  and  if  this  unit  is  interfaced  with  the  HP-41  calculator  they 
reduce  to  1.0/1.0/2.0  whereas  no  data  were  available  in  case  the  printer 
and  the  cassette  drive  are  within  the  loop.  Their  suggestion  with  respect  to 
the  very  slow  data  acquisition  rate  was  to  definitely  take  out  of  the  loop 
the  printer  and  the  cassette  drive  and  also  try  to  use  AUTOZERO  OFF  and 
AUTORANGING,  or  switch  to  another  system  as  the  HP-85  desk  top  computer 
or  the  neuer  HP-71  which  provide  the  opportunity  of  using  the  advanced 
commands  that  allow  for  better  programming  flexibility  (i.e.  use  of 
AUTOZERO  and  AUTORANGING  functions  and  as  less  digits  of  resolution  as 
possible  according  to  the  requirements  of  the  experiment)  and  also  by-pass 
the  problem  of  storing  and  printing  data  by  using  their  own  modules 
without  contributing  to  a  larger  interface  loop.  Under  the  above 
circumstances  (interfacing  the  HP-3421A  DACU  with  the  HP-85  and  using 
the  advanced  commands)  according  to  the  HEWLETT-PACKARD  engineers  the 
recording  rates  of  the  HP-3421A  DACU  as  presented  in  the  Operating, 
Programming  and  Configuration   Manual  [Ref.  15  :  p.  32]  reduce  to    3/9/1  1 
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respectively  depending  on  the  number  of  digits  of  resolution  required  which 
is  substantially  higher  to  what  already  has  been  achieved. 

The  points  of  contact  with  HEWLETT-PACKARD  as  well  as  the 
official  correspondence  made  with  the  customers  division,  by  which  the 
company  was  asked  to  confirm  by  letter  the  information  released  by  phone 
and  suggest  possible  solutions  of  improvement  is  presented  in  Appendix  B. 


C.         OVERALL  EXPERIMENTAL  SET-UP 

To  come  up  with  the  overall  experimental  set-up  several  caclulations 
had  to  be  made  concerning  the  optimum,  yet  safe  number  of  pieces  of 
equipment  that  ought  to  be  used.  These  calculations  are  presented  in 
Appendix  C  and  dictated  the  utilization  of  one  (1)  power  supply  HP-6216B 
for  every  set  of  nine  (9)  Load  Cells  (  per  INSTRON  1000  )  as  well  as  an  extra 
power  supply  HP-6216B  for  the  LVDTs  that  would  be  externally  mounted  to 
the  INSTRON  machines  for  the  recording  of  the  creep-rupture  displacement. 

The  basic  steps  that  were  followed  for  the  overall  experimental 
set-up  which  is  presented  in  a  diagram  form  by  Figures  10a  and  10b,  are 
also  given  in  Appendix  C  under  the  title  "  Implementation  of  the  Overall 
Experimental  Set-Up ". 
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Figure    10a 


Block  diagram  of  the  overall  experimental 
set-up  involving  four  INSTRON  1000  machines 
proposed  for  the  Loading  Sequence  phase. 
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Figure  10b  :  Block  diagram  of  the  overall  experimental 
set-up  involving  four  INSTRON  1000  machines 
proposed  for  the  Creep-Rupture  phase. 
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V.   DI5CU551QN-RE5ULT5 


A.         IN5TR0N  LIMITATIONS 

According  to  the  INSTRON  manual  [  Ref.  14  ]  the  Liquid  Crystal  Display 
(LCD)  indicator  for  the  extension  readout  of  the  moving  cross-head  of  the 
INSTRON  between  the  no-load  case  and  the  under  load  situation,  gives  an 
agreement  within  about  10  %  and  hence  only  the  latter,  which  is  directly 
related  to  the  ultimate  experiment,  was  thoroughly  examined. 

Several  attempts  were  made  to  obtain  displacements  for  the 
minimum  knob  speed  setting  (i.e  0.5  in/min-MIN)  but  it  turned  out  that  no 
movement  of  the  cross-head  had  taken  place  as  it  was  observed  by  both  the 
LCD  and  the  DIAL  indicators;  the  latter  was  used  to  monitor  the  physical 
displacement  as  described  in  the  Appendix  A.  Furthermore  it  was  observed 
that  after  a  short  period  of  time  the  LCD  displacement  readout  started  to 
decrease  with  time. 

Using  the  bracketing  technique  it  was  found  that  no  matter  what  the 
load  was  (within  the  capability  of  the  equipment)  no  stall  occured  for  knob 
speed  settings  larger  than    0.5  in/min-30  %  and  the  overtravel-lim iters 
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system  had  finally  to  interrupt  the  loading  sequence  to  protect  the 
equipment  from  overload. 

A  good  number  of  attempts  were  also  performed  at  the  range  of 
speeds  (0.5  in/min-MIN,  0.5  in/min-12  %).  This  showed  completely 
irrelevant  indications  (Table  I)  as  far  as  the  observed  and  the  expected 
speeds  correlation  was  concerned  and  suggest  that  in  actual lity  the 
cross-head  had  moved  so  little  and  in  such  a  random  way  that  no 
significance  could  be  given  to  these  small  Variable  Speed  Indicator  (VSI) 
settings. 

Several  knob  speed  settings  were  tried  to  evaluate  the  quality  of  the 
response  and  some  representative  sample  data  are  given  in  Table  II  which 
actually  turned  out  to  be  in  very  good  agreement  with  the  anticipated 
results.    One    can    also    observe    from    Figures     11a,    lib    and    lie 

corresponding  to  the  above  mentioned  table  that  the  functions  d  D(A|_  =  f(t) 

d  |_£Q  =  f(t)  and  P  =  f(t)  are  single  valued,  monotonically  increasing  and  as 

the  time  increases  there  is  a  tendency  for  the  displacement  and  the  load 
intervals  to  decrease.  Having  thus  verified  the  anticipated  results  and 
observing  from  the  quality  of  the  response,  that  the  helical  spring  employed 
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TABLE    i 

INSTRON  1000  RESPONSE  FOR  LOW  SETTINGS 
OF  THE  VARIABLE  SPEED  INDICATOR 
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TABLE    II 

TYPICAL  DATA  BASE  SHOWING  THE  QUALITY  OF  THE  RESPONSE 
OF  THE  INSTRON  1000  TESTING  EQUIPMENT 
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Typical  response  for  the  Loading-Sequence. 
Variation  of  physical  displacement  (DIAL) 
of  the  INSTRON  1000  moving  cross  -  head 
with  time. 
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Figure  lib:  Typical  response  for  the  Loading-Sequence. 
Variation  of  displayed  displacement  (LCD) 
of  the  INSTRON  1000  moving  cross  -  head 
with  time. 
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managed  to  represent  quite  well  the  linear  model  assumed,  the  data  base  of 
Table  III  was  created  by  repetition  of  the  loading  sequence  at  various 
speed  settings.  A  closer  look  to  this  data  base  yielded  the  following  results. 
There  is  always  an  one  to  one  correspondence  between  the  following 
pairs  of  variables  of  interest: 


1.  The  LCD   and  the  DIAL  displacements.  Their  relationship  appears 
to  be  a  linear  under  a  slope  of  one  ( 1 )  as  shown  in  Figure  1 2a.  The 

greater  error  (noise)  band  which  is  observed  to  the  lower  side  of 
the  simple  curve  fit  is  attributed  to  the  sequence  under  which  the 
data  were  manually  recorded.  Had  the  data  been  recorded  in  the 
opposite  sequence,  a  greater  error  band  to  the  upper  side  of  the 
simple  curve  fit  would  have  heen  observed. 

2.  The  expected  speed  and  the  observed  speed  as  shown  in  Figure 
12b. 

3.  The  stall  load  and  the  Variable  Speed  Indicator  (VSI)  as  shown  in 
Figure  12c. 


The  plot  of  stall  load  versus  the  VSI  as  presented  in  Figure  12c 
enables  the  descission  maker  to  predict  the  appropriate  knob  combinations 
to  reach  a  certain  load  level  without  stalling.  However  this  plot  is 
characterized  by  a  lot  of  scatter  in  the  data  and  that  is  why  it  is  also 
presented  in  a  form  to  compensate  for  the  mean  upper  and  lower  values 
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TABLE    III 

DATA  BASE  OF  THE  STALL  CONDITIONS  DURING  THE  SIMULATION 
FOR  VARIOUS  SETTINGS  OF  THE  VARIABLE  SPEED  INDICATOR 
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TABLE     III 

DATA  BASE  OF  THE  STALL  CONDITIONS  DURING  THE  SIMULATION 
FOR  VARIOUS  SETTINGS  OF  THE  VARIABLE  SPEED  INDICATOR 
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TABLE    III 

DATA  BASE  OF  THE  STALL  CONDITIONS  DURING  THE  SIMULATION 
FOR  VARIOUS  SETTINGS  OF  THE  VARIABLE  SPEED  INDICATOR 
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values  achieved  at  every  VSI  knob  setting  as  Figure  12c- 1  after  the 
appropriate  data  manipulation  presented  in  Table  lll-A.  It  is  therefore 
evident  that  in  order  to  avoid  stall  one  has  to  take  into  account  the  V5I 
knob  setting  corresponding  to  the  mean  which  also  gives  the  higher 
confidence  level  in  that  no  stall  will  occur.  For  the  actual  experiment 
where  nine  (9)  fiber  bundles  will  have  to  be  tested  simoultaneously  totaling 
450  lbs  order  of  magnitude  load,  Figure  1 2c- 1  suggests  that  the  VSI  knob 
setting  to  be  used  has  to  be  that  of  0.5  in/min-21  %. 

The  graphs  of  Figures  12a,  b,  c,  c-1  are  characterized  by  a  great 
amount  of  scatter  which  appears  to  be  of  increasing  magnitude  as  the  VSI 
setting  is  increased.  However  this  is  not  true  and  this  indication  is  only 
attributed  to  the  limitations  of  the  testing  equipment  itself.  Had  the 
capacity  of  the  1NSTR0N  been  larger  the  scatter  would  decrease  almost  to 
zero.  This  statement  ows  its  credit  to  the  large  amounts  of  attempts  that 
have  been  made  for  VSI  knob  settings  larger  than  0.24  in  which  the  switch 
limiter-system  had  always  taken  over  to  interrupt  the  loading  sequence  for 
safety  purposes  in  the  neighborhood  of  950-975  lbs.  Being  conservative  this 
means  that  a  stall  load  of  at  least  1000  lbs  for  the  VSI  of  0.25  and  1200 
lbs  for  the  VSI  of  0.26  would  have  been  achieved  if  capacity  were  not  a 
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TABLE    lll-A 


MIN,  MEAN  &  MAX  VALUES  OF  PARAMETERS  OF  INTEREST 
FOR  VARIOUS  VSI  SETTINGS  AT  STALL  CONDITIONS 


Column  1 

'  0.13 
0.14 
0.15 
0.16 
0.17 
0.10 
0.19 
020 
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022 
0.23 
024 
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Column  4 
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16.0 

70 

13.1 

19.6 

75 

232 

30.7 

SO 

4.1 

31.5 

85 

372 

433 

90 

40.7 

502 

95 

43.1 

60.6 

100 

373 

33.0 

105 

412 

58.9 

110 

75.0 

802 

115 

69.7 

82.8 

120 

71.4 

90.8 

Column  5 

22.8 
303 
35.4 
42.0 
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TABLE    lll-A 

MIN,  MEAN  &  MAX  VALUES  OF  PARAMETERS  OF  INTEREST 
FOR  VARIOUS  VSI  SETTINGS  AT  STALL  CONDITIONS 


oiumn  6 

Column  7 

Column  8 

Column  9 

Column  10 

Column  1 1 

15.0 

141.9 

638.0 

0.0 

14.5 

65.0 

30.0 

167.3 

411.0 
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17.7 

75.0 

112.0 

714.1 

1529.0 

11.0 
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30.0 
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65.0 
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10.0 
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304.0 

167.0 
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2145.0 

15.0 

176.6 

941.0 

1048.0 
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2207.0 
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TABLE    Ml -A 

MINI,  MEAN  &  MAX  VALUES  OF  PARAMETERS  OF  INTEREST 
FOR  VARIOUS  VSI  SETTINGS  AT  STALL  CONDITIONS 


Column  12 

Column  13 

Column  14 

Column  15 

Column  16 

Column  17 

65.0 

503.4 

1961.0 

62.0 

89.5 

201.0 

141.0 

468 .3 

1885.0 

69.0 
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226.0 

251.0 

14395 
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Response  Quality  of  the  Simulated  Experiment. 
Relations  between  MIN,  MEAN  and  MAX  achieved  stall 
load  and  INSTRON  1000  variable  speed  indicator 
knob    setting    (VSI    KNB    STNG). 
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problem.  Including  these  postulated  values  in  Table  lll-A  one  can  come  up 
with  Figure  12c-2  which  clearly  shows  the  decrease  of  the  error  band 
mentioned  above. 

The  accuracy  of  the  LCD  is  limited  to  O.Ol  in  and  therefore 
displacements  of  lower  order  of  magnitude  cannot  be  correctly  recorded  as 
for  example  during  the  graphite  fiber's  creep.  Roughly  the  deformation 
during  creep  is  about  20%  of  that  achieved  at  the  end  of  the  elastic 
behaviour.  Knowing  that  the  elastic  deformation  is  of  the  order  of  magnitude 
of  2  %,  for  the  gauge  length  of  10  inches  to  be  used,  the  ultimate  creep 
displacement  would  be  of  the  order  of  magnitude  of  0.04  in  and  therefore  it 
is  unlikely  that  the  IN5TR0N  LCD  could  record  any  accurate  data  during  the 
creep  experiment. 

Another  point  worthy  of  comment  is  that  from  the  simulation  of  the 
actual  experiment  via  the  helical  spring,  it  was  observed  that  whenever  a 
stall  load  level  was  reached  the  LCD  display  appeared  to  decrease  in  time, 
the  order  of  magnitude  of  which  was  estimated  to  be  about  0.01  in  /  20  sec. 
Hence  the  LCD  display  cannot  be  considered  a  reliable  means  of  recording 
the  fiber*s  creep  data  which  is  of  primary  importance  for  life  testing.  The 
reason  for  this  decrease  in  the  LCD  displacement  reading  lies  on  the  fact 
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that  this  reading  is  achieved  through  the  Shaft  Encoder  (SE),  a  device  that 
generates  an  output  in  proportion  to  the  revolutions  per  minute  of  the 
Intermediate  Drive  Shaft  (IDS)  by  which  it  is  driven.  Consequently  in  case 
of  a  stall,  due  to  the  internal  friction,  the  IDS  is  not  rotating,  no  output  is 
generated  by  the  SE  and  hence  the  decrease  in  the  LCD  displacement 
reading. 

Finally  summarizing  the  validity  of  these  results  one  can  rely  on  the 
INSTRON  1000  only  for  the  loading  sequence  until  the  stall  load  is  achieved 
by  merely  selecting  VSI  knob  settings  as  dictated  by  the  relatioship  of 
Figure  12c- 1,  but  cannot  rely  on  the  LCD  for  the  creep-rupture  phase  and 
therefore  other  means  of  recording  and/or  displaying  the  displacement 
during  the  creep  experiment  have  to  be  implemented.  Such  means  is  the 
SCHAEVITZ  LVDT  (for  recording)  mounted  on  the  loading  frame  of  the 
INSTRON  in  conjunction  with  a  high  accuracy  multimeter  as  the  available 
FLUKE  8440A  to  measure  and  display  it's  output. 


B.         LVDT  LIMITATIONS 

The  bracketing  technique  outlined  in  Section  III  of  the  Appendix  A 
under  the  title  "Procedure"  showed  that  the  LVDT  was  totally  insensitive 
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for  an  input  of  less  than  1 1.0635  V.  This  result  means  that  for  all  practical 
purposes  the  output  is  too  low  for  being  able  to  be  recorded  and/or 
displayed  for  an  input  less  than  the  12  V  order  of  magnitude.  On  this  basis 
the  LVDT  response  was  then  checked  at  three  points  namely  24  V,  18  V  and 
12  V  respectively,  for  the  displacements  of  0.02  in,  0.04  in  and  0.06  inches 
due  to  the  fact  that  during  creep,  the  order  of  magnitude  of  the  creep 
displacement  for  the  graphite  fiber  bundles  is  known  to  be  of  the  0.04  order 
of  magnitude  (  2  %  of  20  %  of  1 0  in  ). 

For  a  positive  direction  of  displacement  (  i.e  0.02,  0.04,  0.06  in  )  and 
after  averaging  the  output  of  about  10  readings   per  case,  one  can  come  up 

with  Table    IV-A  and  observe  that  ,  the  closer  the  input  voltage  (Vjn)  to 

the  calibration  value  of  the  24  V,  the  better  the  response  and  vise  versa.  In 
fact  the  normalized  response  of  the  LVDT  at  the  end  points  of  the  voltage 
region  (18V-24V)  showed  almost  exactly  the  same  slope,  as  it  was 
anticipated  and  a  high  goodness  of  fit  as  shown  in  Figures  1 3a  and  1 3b. 

For  a  negative  direction  of  displacement  (  i.e  0.06,  0.04,  0.02  )  one 
can  similarly  come  up  with  Table  IV-B  which  leads  to  exactly  a 
similar  observation  as  one  can  see  by  looking  at  Figures  1 3c  and  1 3d. 
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TABLE    IV- A 

LVDT  OUTPUT  FOR  VARIOUS  INPUTS  AND  POSITIVE 
DIRECTION  OF  DISPLACEMENT 


Column  1 

0.02 
0.04 
006 

Coiumn  2 

-1.76213 
-1.71797 
-1.31844 

Column  3 

-0  14684 
-0.14316 
-0.10987 

Column  4 

-3.52404 
-3.34923 
-3.29234 

Column  5 

-0.19578 
-0.19218 
-0.18291 

Column  6 

-5.57715 
-5.46213 
-5.22860 

Column  7 

-0.23238 

HJ  .22739 
h3  .21786 

OSPtCTNT(in) 
PstvOrctn 

AVR6  VoutiV) 
Vin  *  1 2V 

Vout  /  Vln 
Vin»  12V 

AVRSVout(V) 
Vin  •  18V 

Voul/Vin 
Vin  *  1 8V 

AVR6VouUV) 
Vin  »  24V 

V  out  /  Vin 
Vin  »  24V 

TABLE    IV-B 

LVDT  OUTPUT  FOR  VARIOUS  INPUTS  AND  NEGATIVE 
DIRECTION  OF  DISPLACEMENT 


Column  1 

0.06 
0.04 
0.02 

Column  2 

-1  63297 
-1.71633 
-1. 76616 

Column  3 

-0.13608 
-0.14303 
-0.14718 

Column  4 

-3 .30045 
-3.43410 
-3.53410 

ColumnS 

-0.18336 
-0.19189 
-0.19634 

Coiumn  6 

-5.24423 
-3.46877 
-5.59618 

Column  7 

hD.21851 
-0.22737 
-0  23317 

0SP1.CnNT(in) 
Ngtv  OrrLn 

AVR6  Vout(V) 
Vin*  12V 

Vout  /  Vin 
Vin*  12V 

AVRGVout(V) 
Vin  ■  1 8V 

Vout  /  Vin 
Vin*  18V 

AVT36VouKV) 
Vin  »  24V 

V  out  /  Vln 
Vin  »  24V 
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Figure    13a 


Response  of  the  LVDT  for  various  inputs 
and    directions  of  displacement 
Variation  of  Average  Vout  with  Positive 
displacement. 
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Response  of  the  LVDT  for  various  inputs 
and    directions  of  displacement. 
Variation  of  (  Vout  /  Vin  )  with  Positive 
displacement. 


107 


-2- 


(Vin*  12V) 

a 


y  =  -  I .8384  +  3.3298x    R  =  0.99 


-     "3 

a 
e 
> 


> 
< 


-5- 


3  Column  2 
♦  Column  A 
■     Column  6 


(Vina  iav) 


y-  -3.6632  +  5.84 12x    R  -  0.98 


(Vin  =  24V) 


5.7883  +  8.7988x    R  =  0.99 


0.01 


0.02 


0.03 


0.04 


0.05 


0.06 


Hqtv  Displacam*nt  (in)— Column   I 


Figure    13c 


Response  of  the  LVDT  for  various  inputs 
and    directions  of  displacement. 
Variation  of  Average  Vout  with  Negative 
displacement. 
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Response  of  the  LVDT  for  various  inputs 
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Variation  of  (  Vout  /  Vin  )  with  Negative 
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Furthermore  by  comparing  Figures  13b  and   13d  one  can  see  that  for 

the  voltage  region  (18V-24V)  the  difference  in  slope  of  the  corresponding 
curves  is  only  at  the  third  decimal  point  (  less  than  05%) . 

Consequently  for  this  region  the  LVDT  response  is  almost 
independent  of  the  direction  of  displacement. 

The  stability  of  the  LVDT  in  time  was  checked  by  recording  its' 
output  at  a  step  inteval  of  5  seconds  for  about  1  hour  at  the  input  voltage 
check  points  of  12v,  18V  and  24V  as  shown  in  Table  V.  Plotting  of  the 
results  in  Figures  14a,  14b  and  14c  shows  that  although  the  output  is 
pretty  stable  with  time  prior  to  its*  stabilization  overshoots  for  a  period  of 
10  to  20  seconds  depending  on  the  input  voltage  used  (the  lower  the  input 
voltage  the  longer  the  overshoot)  and  finally  stabilizes  rather  quickly  under 
a  damped  oscillation.  It  is  believed  that  this  is  the  reason  that  the  data 
recorded  during  the  initial  run  of  the  experiment  showed  that  the 
cross-head  moved  up  and  down.  Furthermore  the  reason  that  a  polarity  check 
did  not  reveal  any  faults,  is  that  the  amplitude  of  the  oscillation  is  very 
small  (  0.000 IV  )  compared  to  the  output  of  the  LVDT  which  is  of  the  order 
of  several  volts. 

The  check  of  the  time  dependence  was  made  possible  by  monitoring 
the  output  of  the  LVDT  in  the  neighborhood  of  1  hour,  10  hours  and  15  hours 
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TABLE   V 


TYPICAL  LVDT  TIME  RESPONSE  FOR  VARIOUS  INPUTS 


Column  1 

Column  2 

Column  3 

Column  4 

0 

-1.72075 

-3.47150 

-5.46130 

5 

-1.71997 

-3.44820 

-5.46330 

10 

-1.71879 

-3.45970 

-5.46180 

15 

-1.72059 

-3.45760 

-5.46180 

20 

-1.71967 

-3.45760 

-5.46180 

25 

-1.71824 

-3.45760 

-5.46180 

30 

-1.72016 

-3.45760 

-5.46180 

35 

-1.71930 

-3.45760 

-5.46180 

40 

-1.71788 

-3.45760 

-5.46180 

45 

-1.71135 

-3.45760 

-5.46180 

50 

-1.71099 

-3.45760 

-5.46180 

TIKE  (sec  ) 

Vout  (V)  0 

Vout  (V)  9 

Vout(V)   0 

Vinr  12V 

Vin=  18V 

Vin  =  24V 

11 1 


3 
O 
> 


1.710 


•1.712 


■1.714  - 


-1.716  - 


-1.713  - 


-1.720  - 


Vin  =  12V   2,  Pstv  drctn  of  dsplcmnt 


a    Column  2 


-1.722  4—| — i — i—i — i — ■ — i — '  i  ' — i — i — i — i — i — i — i — i — i — i  i  i 
0    5    10   15   20   25   30   35   40   45   50 


TIME  (  sec  )-Column  1 


Figure    !4a       :     Time  response  of  the  LVDT  for  Vin  -  12V. 
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Figure    14b       :     Time  response  of  the  LVDT  for  Vin  -  18V. 
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Figure    14c       :     Time  response  or  Ihc  LVDT  for  Vin  ■  24V. 


for  the  same  displacements  of  0.02,  0.04  and  0.06  inch  in  both  directions 
and  for  the  same  input  voltage  check  points  (12V,  18V  and  24  V).  The  results 
are  presented  in  Tables  V-A  and  V-B  and  when  plotted  as  in  Figures  1 5a 
and  15b  showed  that  the  minimum  time  dependance  (  less  than  10~4 
relative  errror  )  is  achieved  for  both  directions  of  displacement  at  24V  and 
at  about  0.04  in  displacement.  Although  negligible  (  less  than  2x1 0~^ 
relative  error  )  the  time  dependance  at  18V  was  almost  constant  with  a 
shallow  minimum  at  about  0.03-0.055  in,  whereas  no  minimum  point  was 
identified  in  the  12V  case. 

Consequently  for  the  entire  voltage  region. (12V  -24V)  the  LVDT 
response  shows  negligible  time  dependence. 

Since  the  area  of  interest  as  far  as  the  creep  is  concerned  is  the 
region  between  0.020  in  and  0.025  inches,  drifting  in  time  does  not  seem  to 
be  a  serious  problem  for  any  input  voltage  within  the  above  specified  region; 
however  a  slight  dependence  can  be  observed  in  the  lower  half  of  this 
voltage  region  and  hence  it  is  prudent  to  use  the  24V  input  voltage  for  the 
LVDTs  that  are  going  to  be  used. 
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TABLE   V-A 

LVDT  RELATIVE  ERROR  FOR  VARIOUS  INPUTS  AND  POSITIVE 
DIRECTION  OF  DISPLACEMENT 

Column  1  Column  2  Column  3  Column  4 


0.02 
0.04 
0.06 


DSPLCMNT(in) 
Pstv  Drctn 


0.00003© 
0.0008 10 
0.00132© 


DVout  /  Vln 
Vin  =  12V 


0.00128© 
0.00129© 
0.00141© 


DVout  /  Vin 
Vin  =  18V 


0.00154© 
0.00002* 
0.00159* 


DVout  /  Vin 
Vin  =  24V 


CCCES  : 


Ihr 


:10hrs 


TABLE  V-B 

LVDT  RELATIVE  ERROR  FOR  VARIOUS  INPUTS  AND  NEGATIVE 
DIRECTION  OF  DISPLACEMENT 


Column  1 

0.06 
0.04 
0.02 

Column  2 

0.00029* 
0.00051© 
0.00003© 

Column  3 

0.00105© 
0.00111* 
0.00146© 

Column  4 

0.00164© 
0.00007© 
0.00157© 

DSPlCMNT(in) 
Ngtv  Drctn 

DVout  /  Vin 
Vin  =  12V 

DVout  /  Vin 
Vin=  18V 

DVout  /  Vin 
Vin  =  24V 

COOES    : 

*  :15  hrs 

•  :10hrs 

©   :    Ihr 
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Figure    15a 


:     LVDT  relative  error  for  various  inputs  and  Positive 
direction  of  displacement. 
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Figure    15b 


LVDT  relative  error  for  various  inputs  and  Negative 
direction  of  displacement. 


118 


V.   CONCLUSIONS-  RECOMMENDATIONS 


In  identifying  the  characterization  methodology  for  efficient 
strength  to  life  data  generation  by  testing  fiber  bundles  (each  may  consist 
of  ten  thousand  (104)  filaments)  it  was  made  clear  that  in  order  for  one  to 
arrive  to  meaningful  1  statistical  data,  long  periods  of  time  and  a  very 
efficient  data  acquisition  system  is  required.  The  high  efficiency  of  the 
data  acquisition  system  is  primarily  needed  during  the  loading  phase 
because  the  life  prediction  is  based  on  the  conditional  probability  that 
no  failures  have  occurred  during  the  loading  sequence  and  this  probability 
can  be  calculated  only  if  the  number  of  failures  that  occured  during  the 
loading  sequence  is  known;  in  the  stress  rupture  phase,  the  more  accurately 
the  number  of  failures  is  recorded  the  better  the  results.  Furthermore  the 
data  accuracy  during  the  loading  sequence  is  what  gives  rise  to  an  accurate 
lower  tail  of  the  Weibull  life  distribution  which  is  cruisial  in  safe  life 
prediction  calculations. 

The  data  accuracy  during  the  loading  sequence  are  influenced  by  the 
limitations  of  the  testing  equipment  and  the  data  acquisition  system. 
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The  limitations  of  the  testing  equipment  were  examined  by  a 
mechanical  spring  simulation  and  it  was  found  that  in  order  to  avoid 
irregular  indications  during  the  loading  sequence  caused  by  the  stall  of  the 
equipment  due  to  internal  friction,  the  data  base  of  Table  lll-A  and/or 
Figure  12c- 1  have  to  be  used  whereas  for  the  450  lbs  load  level  of 
interest  (simoultaneous  testing  of  nine  (9)  graphite  bundles  Hercules  AS4 
including  the  control  bundle)  the  cross-head  speed  setting  has  to  be  higher 
than    0.5  in  /min  level  with  the  variable  speed  indicator  at  21%. 

The  limitations  of  the  testing  equipment  during  the  creep-rupture 
phase  were  also  examined  by  simulation  and  it  was  found  that  indeed  the 
LVDT  output  is  very  small.  More  specifically,  for  the  recommended 
excitation  voltage  of  the  24V  where  the  time  dependence  errors  are  almost 
inexistent  a  creep  displacement  of  0.  0001  inch  corresponds  to  a  difference 
in  output  slightly  less  than  0.00 IV  which  dictates  that  the  five  digits 
resolution  mode  has  to  be  used  by  the  data  acquisition  system.  The  reason 
that  the  data  recorded  during  the  initial  run  of  the  experiment  was 
considered  unreliable  (in  that  it  showed  that  the  cross-head  moved  up  and 
down,  whereas  only  upward  displacement  is  theoretically  expected  from 
monotonic  creep)  lies  on  the  fact  that  even  though  the    LVDT  output  is 
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stable  with  time,  prior  to  the  composite  sample  stabilization,  it 
overshoots  and  stabilizes  rather  quickly  (within  a  few  minutes)  under  a 
damped  oscillation.  This  oscillation  is  relatively  small  (of  the  0.000 IV 
order  of  magnitude)  compared  to  the  LVDT  output  which  is  of  the  order  of 
several  volts. 

The  data  acquisition  system  used  is  not  capable  of  recording  data  at  a 
rate  faster  than  0.04  readings/ sec.  According  to  the  available  literature 
this  result  appears  to  be  within  the  design  constraints  of  the  instruments 
since  an  interface  loop  involving  only  the  HP-41CX  calculator  and  the 
HP-3421A  DACU  can  barely  give  I  reading/ sec.  It  is  therefore  obvious  why 
the  recording  rate  decreases  substantially  when  the  cassette  drive  and  the 
thermal  printer  are  also  included  in  the  loop,  given  that  they  are  considered 
slow  units  in  the  first  place.  On  the  other  hand  the  data  logger  routine  used 
had  been  modified  to  record  data  at  hundredths  of  seconds  that  is  50  times 
faster  than  the  theoretical  capability  of  the  data  acquisition  unit  itself. 
Finally  due  to  the  limitations  of  the  HP-41CX  calculator  the  advanced 
commands  could  not  be  used  (AUTOZERO  OFF,  AUTORANGING,  etc.)  which  can 
also  substantially  increase  the  recording  rate.  Under  the  above 
circumstances  it  is  recommended  to  substitute  the  existing  data  acquisition 
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system  with  a  simpler  one  composed  of  the  the  HP-342IA  DACU  and  the 
HP-85  desk  top  computer  at  least  during  the  loading  sequence.  The  obvious 
advantages  of  such  a  system  which  can  give  a  reading  rate  per  second  as 
fast  as  3,  9  or  1 1  (for  5,  4  or  3  digits  of  resolution  repectively)  are: 


1.  No  other  peripherals  that  might  limit  the  recording  rates  will 
have  to  be  connected  since  the  HP-85  has  its  own  printer  and 
cassette  drive. 

2.  More  programming  flexibility  will  be  available  through  the  use  of 
the  advanced  commands. 

3.  Substantial  difference  in  memory  capacity  (3 1 9  vs  1 6K). 

4.  During  the  loading  sequence  there  is  no  need  to  record  the  LVDT 
output  due  to  the  relative  rapid  changes,  a  fact  which  saves  one 
channel  and  at  the  same  time  provides  the  opportunity  to  use  the 
three  digits  of  resolution  mode  (N3)  which  can  more  than  triple 
the  recording  rate. 
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APPENDIX    A 

N5TR0N    1000  DESCRIPTION-OPERATING  PRINCIPLES 


A.         GENERAL 

The  INSTRON  model    1000  Universal  Testing  Instrument     shown  in 

Figure  A-l  is  set  to  accept  a  main  power  input  of  120+10  VAC,  single 
phase,  50-60  Hz  and  by  means  of  it's  own  power  transformer,  the  necessary 
AC  voltages  are  provided  to  the  various  systems  of  the  equipment. 

All  the  major  units  of  the  instrument  as  well  as  their  sub-units  are 
described  with  respect  to  Figure  A- 1  by  an  item  number  given  in  parenthesis 
next  to  the  unit  under  description. 


B.         MAJOR  UNITS 

These  can  be  thought  of  being  the  loading  frame  (1),  the  drive  train 
assembly  and  the  control  console  (19). 

1.        Loading  Frame 

This  is     the    part      of     the     equipment     where     the     test 
specimens  are  mounted  and  the  loading  (tension  or  compression)  is  applied. 
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Figure   A-l  :     Major  Units  and  sub  units  of  the  testing  Equipment. 
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2.        Drive  Train 

This    is   a    series   of    several    units    the   prime    function     of 
each  one  can  be  described  as  follows: 


a.  Drive  Motor  (20) 

It's  rotary  motion  is  transfered  to  an  intermediate  drive  shaft  by 
means  of  a  motor  belt  (2 1 ). 

b.  Intermediate  Shaft  (22*) 

Its  rotary  motion  is  in  turn  transfered  to  the  leadscrew  drive 
pulleys  (25)  by  means  of  the  main  drive  belt  (23). 

c.  Leadscrew  Drive  Pulleys  (25) 

These  are  rotated  by  the  main  drive  belt  (23)  and  are  nechanically 
connected  to  the  leadscrews  (6). 

d.  Leadscrews  (6) 

These  are  threaded  rods  which  are  at  the  top  supported  by  the 
fixed  cross-head  structure  (7)  and  can  turn  the  moving  cross-head 
(2)  at  commanded  speeds. 

e.  Moving  Cross-Head  (2) 

This  is  the  element  within  the  loading  frame  ( 1 ),  that  applies  the 
loads  to  the  test  specimen  at  commanded  speeds. 


3.        Control  Console 

This  part  of  the  equipment  serves  as  a  housing  for  the 
controls  and  the  indicator  systems  (19).  A  detailed  view  of  it's  front  and 
rear  panel  is  presented  in  Figure  A-2. 
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Figure   A-2  :     Front  and  rear  Control  Panels  of  the  testing  Equipment. 
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C.         OTHER  IMPORTANT  FEATURES 

Among  the  important  features  of  the  IN5TRON  1000  which  are 
described  in  full  detail  in  the  manual  of  the  equipment  [Ref.  14],  the 
following  two  need  to  be  generally  described: 

1.        Overtravel-Limiters  System 

The  purpose  of  this  system  is  to  protect  the  Load  Transducer 
(LT,  3)  from  overload  by  stopping  the  moving  cross-head  (2).  This  is 
achieved  in  the  following  sequence: 


a.  Setting  the  adjustable  upper  and  lower  limit  stops  (13,14)  in 
appropriate  positions  on  the  limit  switch  rod  (12). 

b.  Activating  the  overtravel  limit  switches  (11)  after  the  limit 
switch  actuator  (15)  which  is  mounted  on  the  moving  cross-head, 
gets  in  touch  with  either  one  of  the  adjustable  stops. 


2.        IMC  Digital  Processor/Controller  model  700  (IMC  DP/C) 

The  purpose  of  this  feature  shown  in  Figure  A- 3  is  associated 
with  the  upgrading  of  the  INSTRON  1000  so  that  testing  of  viscoelastic 
materials  could  also  be  possible.  More  specifically  this  unit  is  capable  to 
create  a  constant  stress  level  in  time,  using  the  applied  load  as  the  control 
variable  which  is  of  prime  importance  during  the  creep.  This  is  achieved  by 
the  following  sequence: 
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a.  Setting  the  desired  maximum  and  minimum  load  levels  on  the  face 
of  the  IMC  DP/C  defining  thus  the  permissible  band  of  the  load 
variation. 

b.  Activating  the  automatic  mode  by  selecting  the  appropriate  button 
on  the  face  of  the  IMC  DP/C  which  allows  the  moving  cross-head 
to  move  under  its  control  to  the  desired  maximum  setting  and  then 
continually  adjust  to  correct  for  the  effects  of  the  creep,  within 
the  range  setting  of  the  previous  step. 


D.         FUNCTIONS 

I.        Load  Measurement 

The  load  measurement  is  achieved  by  means  of  a  calibrated 
Load  Transducer  (  LT,  3  )  that  has  to  be  installed  into  the  moving 
cross-head.  There  are  three  options  in  calibrated  LTs  (10,  100  and  1000 
lbs)  dependent  on  the  anticipated  load  level.  The  input  of  the  transducer  is 
transformed  to  an  electrical  signal  by  means  of  a  LINEAR  VARIABLE 
DIFFERENTIAL  TRANSFORMER  (LVDT).  This  is  a  sensitive  device  that  can 
measure  minute  deflections  of  the  load  transducer  in  the  load  weighting 
system  and  is  composed  of  two  major  elements  shown  in  Figure  A-4: 

a.  The  LVDT  coil  assembly  (2*)  and 

b.  The  LVDT  core  assembly  (4*). 
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Figure    A-4  :     LVDT  coil  and  core  assemblies. 
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It's  transformed  electrical  signal  is  finally  digitally  displayed 
to  a  Liquid  Crystal  Display  (LCD)  panel  located  at  the  upper  right  corner  of 
the  front  side  of  the  console,  at  a  frequency  dependent  on  the  commanded 
cross  head  speed. 

2.  Extension  Measurement 

The  extension  measurement  of  the  moving  cross-head  is 
achieved  by  utilizing  a  SHAFT  ENCODER  (SE,  22),  a  device  which  generates 
an  output  in  proportion  to  the  revolutions  per  minute  of  the  intermediate 
shaft  (  IS,  22*  )  by  which  it  is  driven.  This  output  is  an  incremental 
measure  of  the  displacement  of  the  moving  cross-head  and  it  is  also  used  as 
a  feedback  to  maintain  the  commanded  speed.  The  defficiency  of 
the  extension  measurement  system  lies  on  the  fact  that  no  output  can  be 
generated  by  the  SE  in  any  case  case  the  equipment  could  not  overcome  its 
own  internal  friction  as  in  a  case  of  a  partial  or  total  stall. 

3.  Creep  Monitoring 

This  function  is  achieved  via  the  enhancement  feature  of  the 
IMC  DP/C  which  is  capable  of  creating  a  constant  stress  level  in  time 
using  the  applied  load  as  the  control  variable.  Thus,  by  setting  the  desired 
band  of  load  variation  the  moving  cross  head  moves  under  the  IMC  DP/C 
control  by  continually  adjusting  to  the  effects  of  the  creep. 
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I.   IN5TRQN  1000  LIMITATIONS 


A.        OBJECTIVES 


The  objectives  with  respect  to  the  INSTRON  1000  limitations  were: 


1.  To  evaluate  the  minimum  speed  that  the  INSTRON  can  be  operated 
without  stalling. 

2.  To  convert  this  speed  to  a  combination  of  the  INSTRON  labeled 
controls. 

3.  To  check  whether  this  speed  is  load  dependent  or  in  otherwords 
given  a  physically  verified  speed  check  whether  one  is  able  to 
predict  the  corresponding  stall  load. 

4.  To  find  out  any  other  limitations  that  will  become  apparent  during 
the  simulation  experiment. 


B.         THEORETICAL  MODEL 


A  practical  way  to  present  a  relationship  between  a  load  level  (p)  and 
the  equipment  cross-head  speed  (v),  is  to  follow  the  simple  relations 
presented  below  where  t  represents  time  and  1  a  physically  measured 
displacement: 
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(A-1)  (p/t)*(  1/p)  =  1/t 

(A-2)  1/t  =  v 

(A-3)  p  =  v*[  1/(l/p)]*t 


Given  that  the  ratio  (1/p)  can  be  thought  of  representing  the  elastic 
constant  of  a  linear  model  one  could  anticipate  the  following  results: 


1.  The  functions  of    1  =  f(t)  and  p  =  f(t)  to  be  single  valued    and 
monotonically  increasing. 

2.  As  load  increases  the  displacement  and  time  intervals  have  to 
show  a  tendency  to  decrease. 


C.         MECHANICAL  SET-UP 

A  helical  spring  as  the  one  shown  in  Figure  A-5  was  decided  to  be 
used  for  the  simulation  experiment  provided  that  both  the  cross-head  speed 
(v)  and  the  time  (t)  could  be  easily  recorded,  because  it  was  cosidered  the 
best  solution  from  the  safety  standpoint. 

Under  the  requirements  of  comforming  with  the  necessary  safety 
precautions  for  the  operator  and  the  equipment  limitations,  the  mechanical 
set-up  shown  in  Figures  A-6   and   A  7  was  designed  which  provided  the 
capability  of  easily  achieving  the  target  load  level  of  the  450  lbs  order  of 
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HELICAL      SPRING 


DIMENSIONS  :   mates 
SCALE  :      'A 


,        _     &0.36S 


Figure   A-5  :     Helical  Spring  employed  for  the  Simulation. 
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magnitude.  On  the  other  hand  since  one  could  not  rely  on  the  INSTRON 
cross-head  control  knob  markings  for  the  evaluation  of  (v)  due  to  inherent 
limitations  of  the  equipment  [Ref.  14]  and  therefore  a  dial  indicator  (DIAL) 
capable  to  measure  displacement  in  0.001  of  an  inch  (  ten  times  better 
resolution  than  that  of  the  LCD  meter  )  was  also  used  during  the  simulation 
procedure  as  shown  in  Figure  A-8. 

The  overall  mechanical  set-up  for  the  simulation  of  the  loading 
sequence  is  shown  in  Figure   A-9. 


METHODOLOGY 

The  methodology  followed  consisted  basically  of: 


1.  Evaluating  the  quality  of  the  response  as  outlined  in  steps  B.  1  and 
B.2  above. 

2.  Determining  the  lower  threshold  speed  that  the  INSTRON  could  be 
operated  without  producing  irrelevant  indications  and  the  stall 
load  corresponding  to  each  one  of  the  selected  speeds,  using 
the  method  of  bracketing  between  selected  upper  and  lower 
bounds. 

3.  Observing  the  quality  of  agreement  between  the  DIAL  indicator 
utilized  and  the  LCD  meter. 
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Figure    A-8 


Dial  indicator  employed  to  measure  the  physical 
displacement  of  the  IN5TR0N  1000  moving 
cross-head. 
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Figure    A-9  :     Overall  Mechanical  Set-Up  during  the  Simulation. 
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E.         DETAILED  PROCEDURES 


The  detailed  procedures  used  to  achieve  the  data  base  available  are 
presented  below  in  a  check  list  form  to  provide  for  ease  in  reproducing 
parts  or  the  totallity  of  the  results,  under  the  desired  incremental  values 
depending  on  the  application  of  interest. 

1.        Familiarization 

Make  sure  that  you  are  thoroughly  familiar  with  the  controls  of 
the  INSTRON  as  well  as  that  you  have  already  read  and  comprehended  the 
following  istructions. 


CAUTION 

DO  NOT  ATTEMPT  TO  OPERATE  THE  INSTRON  UNLESS 
YOU  ARE  FEELING  COMFORTABLY  WITH   ITS  CONTROLS 
AND  THE  FOLLOWING  INSTRUCTIONS 


No-load  case 


a.  Identify  the  Load  Transducer  on  the  cross-head. 

b.  Calibrate  the  INSTRON  accordingly  [Ref.  1 4 :  ch.  5] 

(1)      Go  through  steps  1   and  2  of  the  calibration  procedure 
(page  5-2  ). 
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(2)  If  using  the  100  lb  or  the  10  lb  capacity  Load 
Transducers  go  through  steps  3  ,4  and  5  (pages  5-2  & 
5-3). 

(a)  For  the  100  lb  capacity  Load  Transducer  go 
through  step  6a  (page  5-4)  anyway,  and  NOTE 
(page  5-4)  as  well  as  steps  7b,  7c  (page  5-5) 
only  if  necessary. 

(b)  For  the  10  lb  capacity  Load  Transducer  go 
through  step  6b  (page  5-4). 

(3)  If  using  the  1000  lb  capacity  load  transducer  go  through 
NOTE  and  step  6a  (page  5-3)  anyway  ,  and  step  7a  (page 
5-4)  only  if  necessary. 

c.  Set  cross-head  control  knob  to  10/.5. 

d.  Set  variable  speed  knob  to  0.1  (first  position  next  to  MIN)  and 
observe  illumination  of  the  corresponding  red  light. 

e.  Set  physical  displacement  scale  (DIAL)  to  starting  reference. 

f.  Based  on  the  settings  of  the  previous  two  steps  make  a  rough 
calculation  of  what  the  maximum  displacement  will  be  according 
to  the  formula: 

(  DISPLACEMENT  )  =  (  SPEED  )  X  (  TIME  ) 

and  set  the  upper  limit  stop.  TIME  has  to  be  assumed  that  of  the 
step  2j  for  the  No-Load  case  and  3c  for  the  load  conditions.  The 
lower  limit  stop  has  to  be  set  so  that  it  matches  with  the  physical 
displacement  starting  reference. 

g.  Set  stop  watch  to  starting  reference  . 
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m. 

n. 
o. 


Make  sure  that  the  appropriate  record  equipment  (lists,  pencil, 
eraser)  is  handy  in  front  of  you  and  that  you  have  understood  how 
the  data  recording  will  be  performed. 

Simultaneously  start  the  watch  timer  and  push  the  white  UP 
button  (tension)  and  observe  it's  illumination  as  well  as,  that  red 
STOP  button's  light  is  turned  off. 

Record  physical  and  LCD  displacements  at  one  (1)  minute 
intervals  for  the  first  ten  (10)  minutes,  at  ten  (10)  minute 
intervals  for  the  first  hour;  if  the  INSTRON  will  stop  earlier  do  not 
forget  to  record  stall  time  and  load  level. 

Push  red  STOP  button  and  observe  it's  illumination  as  well  as  that 
the  white  UP  buttons'  light  is  turned  off. 

Push  yellow  RETURN  button  and  observe  it's  illumination  as  well 
as  that  the  red  STOP  buttons'  light  is  turned  off. 

CAUTION 

THE  RETURN  SPEED  IS  ALWAYS  20  in/min  (500  mm/min) 
AND  THEREFORE  BE  ALLERT  TO  PUSH  THE  STOP  BUTTON 
IN  CASE  ANY  FAILURE  OCCURS  TO  THE  LOWER  LIMIT  STOP 

Set  variable  speed  knob  to  0.2  (  second  position  next  to  MIN  )  and 
observe  illumination  of  the  corresponding  red  light. 

Repeat  steps  2e  up  to  21  above. 

Set  the  variable  speed  knob  to  0.3  (  third  position  next  to  MIN  ) 
observe  the  illumination  of  the  corresponding  red  light  and  repeat 
step  2n  above.  If  the  INSTRON  will  not  stall  do  not  exceed  2/3  of 
the  Load  Transducer  capacity. 

Obtain  plots  of  the  physical  displacement  versus  time  with 
parameter  the  four  different  speed  settings  outlined  in  steps  2c,  d 
m  and  o. 
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3.        Under  load  case 

a.  Repeat  steps  a  to  i  of  paragraph  2.  (No-Load  case). 

b.  Repeat  step  f  of  paragraph  2. 

c.  Record  DIAL  and  LCD  displacements  as  well  as  force  at  one  half 
(1/2)  minute  intervals  for  the  first  ten  (10)  minutes;  hopefully 
the  IN5TRON  will  stall;  do  not  forget  to  record  stall  time  and 
load 

d.  Repeat  steps  2k  and  21. 

e.  Repeat  in  the  order  given  the  following  steps:  2m,  2e,  2f,  2g,  2h, 
2i,  3c,  2k  and  21. 

f.  Repeat  in  the  order  given  the  following  steps:  2o,  2e,  2f,  2g,  2h, 
2i,  3c,  2k  and  21. 

g.  Obtain  plots  of  the  physical  displacement  and  the  applied  load 
versus  time  with  parameter  the  four  diffemt  speed  settings 
outlined  in  steps  2c,  d,  m  and  o. 
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CREEP-DISPLACEMENT/STABILITY  OF  THE  LVDT 


BACKGROUND 


It  is  known  that  for  the  graphite  fiber  bundles  the  order  of  magnitude 
of  the  displacement  during  a  normal  loading  sequence  at  a  constant  stres 
rate  R  is  about  2%  (0.02).  It  is  also  known  that  the  graphite  fiber  bundles 
are  viscoelastic  and  therefore  they  creep  and  the  total  creep  displacement 
is  known  to  be  of  the  0.004  order  of  magnitude  or  in  otherwords  20%  of  the 
maximum  displacement  achieved  at  the  end  of  the  loading  sequence.  Due  to 
the  small  numbers  that  one  has  to  deal  with,  during  the  creep  experiment 
for  a  gauge  length  of  1  =  10  in  it  is  necesary  to  use  a  very  accurate  device 
that  will  be  capable  of  recording  and/or  displaying  the  minute  differences 
in  displacement  during  the  creep.  This  is  the  point  where  the  (LVDT)  comes 
into  play  because  the  minute  creep  displacements  of  the  cross-head  of  the 
IN5TR0N  can  be  converted  to  voltages  of  the  0.1  mV  order  of  magnitude 
which  fortunately  can  be  recorded  and/or  displayed  using  high  accuracy 
multimeters.  And  since  it  is  this  creep-rupture  behavior  that  can  lead  to 
strength-life  data,  the  need  to  be  familiar  with  the  problems  that  may  arise 
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like  the  one  according  to  which  no  output  has  been  recorded  although  the 
IN5TRON  cross-head  had  moved,  is  therefore  obvious.  Some  possible 
explanations  of  the  problem  just  described  could  be  either  the  very  small 
magnitude  of  the  voltage  output  or  drifting  in  time. 


B.  OBJECTIVES 

Therefore  the  objective  is  to  test  the  stability  of  the  LVDT  with 
respect  to  the  following  variables: 

1.  The  Input  Voltage  (Vin). 

2.  The  Direction  of  Displacement  (d-0  or  (d-). 

3.  The  Time  (t) . 

4.  The  Temperature  (T). 

Based  on  the  fact  that  the  actual  experiment  will  be  conducted 
essentially  in  constant  room  temperature  ie.  within  temperature 
differences  of  less  than  five  (5)  degrees,  the  stability  of  the  LVDT  with 
respect  to  the  temperature  does  not  appear  to  be  significant. 


C.         PROCEDURE 

To  implement  the  above  mentioned  objective  the  set-up  presented  in 
Figure  A- 10  was  used  which  is  composed  of  the  following  components: 
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!.  A  regulated  position  test  stand  for  the  mounting  of  both  parts  of 
the  LVDT.  The  core  was  mounted  on  a  micrometer's  vernier  by 
which  a  fully  controlable  insertion  of  the  core  at  increments  of 
O.OOl  in  was  possible. 

2.  A  power  supply  (HP-1616B)  by  which  the  control  of  the  input 
voltage  Vin  was  possible. 

3.  A  high  accuracy  multimeter  (FLUKE  8840A)  for  the  exact 
measurement  of  the  output. 


The  maximum  input  of  the  SCHAEVITZ  LVDT  is  24V  and  therefore  for 
time  saving  purposes  the  "  bracketing  technique  "  was  used  between  the 
limits  of  24V,  12V  and  IV  respectively.  If  this  technique  would  show  that 
some  voltage  region  was  immaterial,  bracketing  would  have  been  then 
applied  for  the  remaining  voltage  region  of  interest. 

For  positive  direction  of  displacement  the  output  voltage  of  at  least 
three  (3)  points  was  obtained,  for  various  selected  values  of  the  input 
voltage  (  at  least  3  ),  and  repetition  of  the  same  procedure  for  the  negative 
direction  of  displacement  was  accomplished. 

The  anticipated  results  according  to  the  specifications  were: 

1.        The  output  of  the  LVDT  to  be  independent  of  Vjn  . 

Check  :  The  ratios  of    Vout/Vin  versus  d,  plot  under  the  same 
slope. 
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The  output  of  the  LVDT  to  be  independent  of  d+  or  d-  . 

Check  :  The  curves  of  Vout/Vjn  versus  d+    and  Vout/Vjn  versus 

d-  coincide. 

The  linearity  of  the  LVDT. 
Check  :  Curves  are  straight  lines. 
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APPENDIX     B 

1.         POINTS  OF  CONTACT 

The  Hewlett-Packard  engineers/technicians  who  served  as  points  of 
contact  in  the  recording  rate  optimization  problem  of  the  HP-3421A  DACU 
and  their  telephone  numbers  are  listed  below  for  further  contacts  if 
necessary: 


1 .  Tom  Turney 

2.  Russ  Mc  Ugh 

3.  Ashley  Werrpass 


(408) 988-7304 
(415)857  8241 
(415)857-5673 


2.         OFFICIAL  CORRESPONDENCE  WITH  HEWLETT-PACKARD 

According  to  HP-IL  Module  OWNER'S  MANUAL  [Ref.  17  :  p.  63]  the 
following  official  letter  was  forwarded  to  the  Customer  Support  Division  of 
the  Hewlett-Packard  Co.  The  purpose  of  this  letter  was  to  explain  the 
existing  situation  with  respect  to  the  slow  data  acquisition  rate  of  the 
HP-3421A  DACU,  have  the  Company  confirm  the  information  released  by 
phone  about  the  recording  rate  capabilities  of  the  equipment  depending  on 
the  peripherals  used  within  the  HP-IL  and  finally  have  it  propose  viable 
solutions  utilizing  the  already  available  equipment. 
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DEPARTMENT  OF  THE  NAVY 

MAVAL  POSTGRAOUATT  SCHOOL 
MONTEREY.  CA  93943-3IOO 


i  *epi_r  °zrzm  TT3 


Curriculum  of  Aeronautical  Engineering 
Major  (HAF)  D.M.  Petndis,  Code  31 
Monterey,  CA  93943-5000 
November  the  21st  1986 


Dear  Sirs, 

I  am  an  officer  of  the  Greek  Air  Force,  currently  stationed  at  the  U.S 
Naval  Postgraduate  School,  for  a  Master  of  Science  Degree  in  Aeronautical 
Engineering.  My  thesis  research  is  related  to  the  graphite  fibers  tensile 
testing,  in  an  effort  to  obtain  a  strength-life  data  base  that  will  serve  in 
predicting  the  service  life  of  a  comoonent  made  of  graphite  fiber 
composites,  given  that  the  strength  of  the  fibers  is  known. 

The  experimental  set-up  is  composed  of  a  tensile  testing  equipment 
(IN5TR0N  1000)  where  the  graphite  fiber  bundles  are  mounted.  Failures  of 
individual  fibers  within  the  bundles  are  transfered  to  an  HP-342lA'0ata 
Acquisition  Control  Unit  by  means  of  a  set  of  high  precision  load  cells 
(INTERFACE  SSM-100).  Installed  in  this  unit  are  two  10  cnannel  multiplexer 
boards  (  44462A-0PT10N  020  ).  During  the  experiment,  data  are  recorded 
using  the  HP-IL  communications  circuit  as  shown  in  the  attached  block 
diagram.  The  HP-44468A  Data  Acquisition  Control  Pac  provided  the  basic 
data  logger  software  program.  This  program  was  modified  to  permit  the 
recording  of  time  in  hundredths  of  seconds  and  edited  to  reduce  the  number 
of  data  registers  used.  A  copy  of  this  program  and  Its  output  are  also 
attached  to  this  letter.  Contrary  to  what  is  advertised  in  page  32  of  the 
Operating,  Programming  and  Configuration  Manual  of  the  3421 A  the  reading 
rates  obtained  are  of  the  0.04  order  o(  magnitude  (see  attacned  plot). 

A  preliminary  telephone  contact  with  your  offices  in  Palo  Alto  on 
September  the  29th  [(415)  857-8175]  provided  us  with  the  information 
according  to  which  the  reading  rates  reduce  to  N5(1)/N4(2)/N3(2)  if  the 
HP-41CX  is  interfaced  to  the  342!  a  as  snown  in  the  attached  block  diagram 
and  can  be  substantially  improved  to    N5(3)/N4(9)/N3(  1 1)  respectively,  in 
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case  the  calculator,  the  thermal  printer  (HP-32I62A)  and  the  digital 
cassete  drive  (HP-82I6IA)  are  substituted  oy  the  HP-85A  desk  top 
computer. 

Similar  information  was  provided  oy  your  offices  in  Santa  Clara 
[(408)  988-7304J  on  OctoDer  the  27th,  and  Palo  Alto  [(415)  857-3241/5673] 
on  OctoDer  the  31st  and  suggestions  to  use  programming  utilizing  the  so 
called  advanced  commands  (AUTOZERO  OFF  and  the  less  digits  of  resolution 
possible)  were  also  made  . 

Based  on  the  information  provided  so  far,  you  are  kindly  requested  for 
the  following: 

1.  To  confirm  by  letter  the  above  information  and  if  possible  provide 
us  with  any  insight  associated  with  the  substantial  reduction  of  the  reading 
rates  in  case  the  3421 A  is  interfaced  with  either  the  hand-held  calculator 
or  the  desk-top  computer. 

2.  To  make  any  helpful  comments  on  the  modified  HP-44468A  Data 
Logger  Routine  and  also  inform  us  whether  any  faster  reading  rates  cauld  be 
obtained  by  reducing  and/or  bypassing  a  number  of  the  HP  devices  in  the 
loob;  if  this  Is  the  case  what  are  your  suggestions  for  the  optimum  scheme? 

3.  To  inform  us  about  the  reasons  that  will  improve  the  recording 
rate  in  case  the  HP-85A  is  interfaced  to  the  3421 A  via  the  HP-IL 

4.  What  increase  of  recording  rate  can  be  achieved  through  the 
advanced  commands  with  either  one,  the  HP-4ICX  calculator  or  the  HP-85A 
computer. 

Your  early  reply  will  be  greatly  appreciated.  In  the  mean  time  we 
remain, 
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Sincerely  yours 


Major  CHAF)  D.M  .Petndis 


PS  In  case  you  have  any  questions  please  call  the  Naval  Postgraduate 

School     Aeronautical     Engineering       Curriculum       Office       at 
(408)  646-2491/2  and  leave  your  message. 

If   later  than     December  the    15th    1986,  please  forward  your 
answer  to  my  thesis  advisor  in  the  following  address  : 

U.S  Naval  Postgraduate  School 
Department  of  Aeronautical  Engineering 

Dr.  Edward  M.  Wu,  Code  67Wt 
Monterey,  CA  93943-5000 
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HP-32162A 
THERMAL  PRINTER 
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APPENDIX  A 


DATA  LOGGER 

01  L3L  "LOSM" 

02  SF  10 

03  GTO  35 

04  LBL  "DLrt" 

05  CF  10 

06  LflL  33 

07  SF  12 

08  SF  27 

09  "  HP  3421 A" 

10  AVIEW 

11  CLKY 

12  "DATA  LOGGER' 

13  AVIEW 

14  cr   12 
is  ••Dt_n" 

16  32 

17  AK 

18  "CLKY" 

19  -32 

20  AK 

21  "LOG" 

23  AK 

24  CF    13 

25  RCL  35 

26  ENTER 

27  INT 
23  X»Y*7 

29  GTO  31 

30  x  i >y 

31  FRO 

32  1  E3 
ZZ-    * 

34  STO  35 

35  LBL  31 

36  "NEW  ?  Y/N" 

37  AON 

33  FC7    10 

39  PROMPT 

40  AOFF 

41  ASTO  Y 

42  "Y" 

43  ASTO  X 


PROGRAH   LISTING 

44  X  =  Y 

45  GTO  01 

46  "EDIT  ?  Y/N" 

47  AON 

43  FC7  10 

49  PROMPT 

50  CF  10 

51  AOFF 

52  ASTO  X 

54  ASTO  Y 

55  X=Y? 

56  GTO  04 

57  SF  12 

53  "— «-E3IT0R— •• 

59  AVIEW 

60  CF  12 

61  GTO  18 

62  LSL  F 

63  "DCV" 

64  AVIEW 

65  1 

66  GTO  02 

67  LSL  02 
63  1  E2 

69  * 

70  1  E-7 

71  * 

72  RCL  05 

73  + 

74  I  So  35 

75  ABS 

76  SF  25 

77  RCL  IND  35 
73  FSTC  25 

79  GTO  10 

30  DSE  35 

31  ABS 

32  "OUT  OF  ROOM1 

33  AVIEW 

34  PSE 

35  F5?  10 

36  GTO  16 


37  GTO  04 

38  LBL  10 

39  RON 

90  FST  10 

91  GTO  21 

92  STO  IND  35 

93  GTO  03 

94  LBL  21 

95  STO  04 

96  RCL  06 

97  39 

98  + 

99  STO  08 

100  1  E-3 

101  * 

102  RCL  35 

103  1 

104  - 

103  STO  07 

106  LASTX 

107  + 
103  * 

109  STO  06 

110  LBL  22 

111  RCL  IND  07 

112  STO  IND  06 

113  DSE  07 

114  ABS 

115  DSE  06 

116  GTO  22 

117  RCL  04 

113  STO  IND  08 

119  RCL  08 

120  INT 

121  38 

123  STO  06 

124  GTO  03 

125  LBL  01 

126  CLX 

127  STO  30' 
123  37 

129  STO  35 


Figure  A.!.    Modified  HP44468A  Data  Logger  Routine 
35 
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130  LB1_  03 

131  FIX  0 

132  CF    29 

133  "FIRST  CH7" 

134  cr   22 

135  AOFF 

136  PROMPT 

137  FC  22 
133  STO  04 

139  ENTER 

140  "LAST  CH?" 

141  PROMPT 

142  CLA 

143  ARCL  Y 

144  ■! —  ' 

145  ARCL  X 

146  AVIEW 

147  1000 
143  / 

149  + 

150  STO  05 

151  CF  23 

152  "FUNCTION?" 

153  PROMPT 

154  FS"  23 
133  GTO  19 

136  "PRESS  FN  KEV 
157  PROMPT 
153  L3L  19 

159  CF  22 

160  "USER0-33  ?" 

161  FRCMFT 

162  FS?  22 

163  GTO  20 

164  "LAST  SETUP" 

165  AVIEW 

166  PSE 

167  "NOT  STORES" 
163  AVIEW 

169  PSE 

170  GTO  03 

171  LBL  20 

172  INT 

173  ABS 

174  34 

175  X<=Y7 

176  GTO  19 

177  RDN 


173  "USER" 
179  ARCL  X 
ISO  AVIEW 

131  16 

132  * 

133  GTO  02 

134  LBL  04 

135  F3?  10 

136  GTO  16 

137  37 

133  RCL  35 
139  X<=Y? 

190  STOP 

191  1  E-3 

192  * 

193  33 

194  + 

195  STO  35 

196  CF  09 

197  "RECORD  Y/N" 

198  CF  23 

199  AON 

200  PROMPT 

201  CF  20 

202  FC723 

203  SF  20 

204  ASTO  Y 

205  "Y" 

206  ASTO  X 

207  X=Y? 
203  SF  09 

209  "RECORD" 

210  FC?  09 

211  "i-OFF" 

212  AVIEW 

213  FC  53 

214  GTO  06 

215  "PRINT?  Y/N" 

216  CF  23 

217  AON 
213  PROMPT 

219  CF  21 

220  FC?  23 

221  SF  21 

222  ASTO  Y 

-^^■P    ■«  w  If 
MW     T 

224  LAST  X 

225  Cf    03 


226  X=«Y? 

227  SF  OS 
223  SF  21 

229  FC?  03 

230  CF  21' 

231  "PRINT" 

232  FC?  03 

233  "i-OFF" 

234  AVIEW 

235  LBL  06 

236  AOFF 

237  SF  29 
233  SF  23 

239  XROM   "ALM" 

240  X<0? 

241  GTO  05 

242  PWRDN 

243  OFF 

244  LSL  05 

243  XEQ  "DLMLM" 

246  ISG  32 

247  GTO  29 

248  1 
*"*4<9  ••  j  •• 

230  SF  25 

251  CLAL 

252  CF  25 

253  PWRDN 

254  OFF 
235  LBL  29 

256  RCL  34 

257  X=0^ 
253  GTO  05 

259  PWRDN 

260  OFF 

261  GTO  05 

262  OFF 

263  LBL  13 

264  SF  10 

265  CF  29 

266  FIX  0 

267  36 
263  RCL  33 

269  X<=Y? 

270  37 

271  STO  35. 

272  LBL  16 

273  "COMMAND  ?" 


Figure  A.  I.    Modified  HP44463A  Data  Logger  Routine  (cont'd) 
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274  AON 

322  RCL  35 

370  GTO  16 

275  FRCMFT 

323  XOY 

371  STO  02 

276  AVIEW 

324  X>Y? 

372  FIX  0 

277  ASTO  X 

325  GTO  13 

373  1  E-3 

278  "LIST" 

326  X=Y^ 

374  * 

279  ASTO  Y 

327  GTO  36 

375  33 

280  X*Y? 

323  XOY 

376  + 

231  GTQ  07 

329  1 

377  STO  35 

282  "INSERT" 

330  - 

373  LBL  08 

233  ASO  Y 

331  1  E-3 

379  RCL  IND  35 

234  x=y> 

332  * 

380  3T0  01 

235  GTO  09 

333  + 

331  XRQM  "DECODE" 

286  "DELETE" 

334  STO  06 

332  ASTO  00 

287  ASTO  Y 

335  LBL  15 

383  CLA 

238  x=y 

336  RCL  06 

334  RCL  35 

239  GTO  12 

337  1 

335  INT 

290  "END" 

333  + 

336  33 

291  ASTO  Y 

339  RCL  IND  X 

337  - 

292  x=y> 

340  STO  IND  06 

338  ARCL  X 

293  GTO  17 

341  ISG  06 

339  "!-:  " 

294  "HELP" 

342  GTO  15 

390  RCL  IND  35 

295  ASTO  Y 

343  LBL  36 

391  ENTER 

296  X=Y"> 

3*4  DSE  35 

392  INT 

297  GTO  14 

345  ABS 

393  ARCL  X 

298  "INVALID  CMD" 

346  GTO  16 

394  RDN 

299  AVIEW 

347  LBL  14 

395  FRC 

300  PSE 

348  "LI3T~ 

396  1  E3 

301  GTO  16 

349  AVIEW 

397  * 

302  L3L  17 

350  PSE 

398  INT 

303  "*END  EDITOR*" 

351  "INSERT" 

399  "! — " 

304  AVIEW 

352  AVIEW 

400  ARCL  X 

305  PSE 

353  PSE 

401  ■•:-,  " 

306  CF  10 

354  "DELETE" 

402  ARCL  00 

307  GTO  04 

355  AVIEW 

403  AVIEW 

308  LBL  12 

356  PSE 

404  ISG  35 

309  CF    22 

357  "END" 

405  GTO  08 

310  "NUMBER  ?" 

358  AVIEW 

406  RCL  02 

311  AOFF 

359  PSE 

407  STO  35 

312  PROMPT 

360  GTO  16 

408  GTO  16 

313  FC7  22 

361  LBL  13 

409  LBL  09 

314  GTO  16 

362  "NONEXISTENT" 

410  CF  22 

315  INT 

363  AVIEW 

411  AOF~ 

316  ABS 

364  PSE 

412  "AFTER  NUnR  "?" 

317  "DELETE" 

365  GTO  16 

413  PROMPT 

318  ARCL  X 

366  LBL  07 

414  FC7  22 

319  AVIEW 

367  37 

415  GTO  16' 

320  33 

368  RCL  35 

416  X<07 

321  + 

369  X<=Y? 

417  -1 

FigureA.1.    Modified  HP44-463A  Data  Logger  Routine  (cont'd) 
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418 

RCL  35 

419 

38 

420 

- 

421 

XOY 

422 

X>Y? 

423 

GTO  13 

424 

STO  06 

425 

"AFTER 

426 

ARCL  X 

427 

AVIEW 

423 

GTO  03 

429 

END 

Figure  A.!.    riodified  HP44468A  Data  Logger  Routine  (con'd) 
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DATA  LCGGER  OUTPUT 

01  L3L  "DLMLrt" 

02  "C" 

0-3  CF  29 
0*  FCT  08 

05  CF  21 

06  RCL  35 

07  STO  01 

08  2 

09  L3L  10 

10  3 

11  + 

12  RCL  IND  01 

13  ENTER 

14  FRC 
13  1  E3 

16  * 

17  INT 

18  XOY 
19 
20 
21 


FCR.1AT 
48 
49 

50 
51 


INT 


22  ISG  01 

23  GTO  10 

24  FIX  0 

25  CLA 

26  ARCL  32 

27  FCT  09 
23  GTO  03 
29  SF  25 

CREATE 
FS-C  25 
GTO  03 
PURGE 
CREATE 
LBL  03 
"PASS  " 
37  ARCL  32 
33  AVISW 

39  CLA 

40  ARCL  32 
0 
FS?  09 

43  3EEKR 

44  CLA 

45  DATE 

46  STO  00 

47  FIX  4 


30 

31 
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41 
42 


34 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

32 

83 

34 

35 

86 

87 

38 

39 

90 

91 

92 

93 

94 


AOATE 
TIME 
STO  01 
FIX  6 
AT  I  ME 
FC?  08 
flVIEW 
PRA 
ADV 
.001 
FS?  09 
WRTRX 

mi 

STO  33 

CF  10 

L3L  02 

RCL  IND  35 

STO  01 

XROM  '•DECaDE'* 

ASTO  00 

ASHF 

ASTO  36 

*" 

4m 

STO  37 

L3L  04 

FIX  0 

CLA 

RCL  01 

XEQ  IND  36 

FS~C  10 

ASTO  00 

STO  IND  37 

CLA 

RCL  01 

INT 

ARCL  X 

RCL  IND  37 

ISG  37 

SIGN 

■•  I  _.     It 

ENG  5 
ARCL  X 

ARCL  00 
FC?  08 
AVIEW 

PRA 


95 

CF  10 

96 

ISG  01 

97 

GTO  04 

98 

CLA 

99 

FIX  0 

100 

ARCL  32 

101 

RCL  ZZ 

102 

FS?  09 

103 

SEEKR 

104 

RCL  IND 

105 

STO  01 

106 

RCL  37 

107 

ST*  33 

108 

1 

109 

- 

110 

1  E3 

111 

/ 

112 

FS?  09 

113 

WRTRX 

114 

ADV 

115 

ISG  33 

116 

GTO  02 

117 

RCL  35 

118 

FRC 

119 

38 

120 

+ 

121 

STO  33 

122 

RTN 

123 

END 

35 


Figure A.Z    notified  HP44463A  Output  Format  Routine 
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AS4  GRAPHfTE,   SPOOL  019 


LEGEND 
a    3UNQLE  -V/SIL'CCME 

*     5L.NCL-  ZCNT^CL 


6  3 

TIME    (MIN) 


o  "°- 
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o  •»- 
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3 3 3 3 3         0         1          J g g 9 n         g       ■  p         3         - 

°- 

>"     ■     — 5 J          -*          -»           '          3 3 3          j          j          3 3 J          j s -j          ^ - 
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Figure  3.2a     Test  2:    Lcad-Time  curves.     Maximum  constant  lead  on  the  central 
sampl8W8s  16.003  kilograms. 
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APPENDIX     C 


I.         ESTIMATION  OF  POWER  SUPPLIES  NEEDED 
a.        Load  Cells  (LC)  Case 

The  Specifications  of  the  HP-6216B  give  a  voltage  range  of 
0-25 V  DC  and  a  current  range  of  0- 400mA,  whereas  the  recommended 
excitation  voltage  of  the  INTERFACE    LCs    is  10V  DC  and  their  resistance 

350  i   3.50. 

A  parallel  solution  evidently  does  not  work  as  it  can  be  easily 
observed  by  Figure  C-1  where  the  current  through  the  loop  exceeds  what 
is  dictated  by  the  power  supply  specifications  for  the  recommended 
excitation  voltage  of  the  LCs  as  shown  below: 


ls  =  vsREC/R ,  2  =  VSREC/  I2j= !  ( 1  /R)  =  1 OV/  2*(  1  /38.888Q) 
ls  =  0.514  A  >  0.400  A 


A  series  solution  involving  two  (2)  INSTRONs  in  which  nine  (9) 
LCs  are  connected  in  parallel  as  shown  in  Figure  C-2  suggests  that  no 
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I 
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J 

i          1 

!             1 

I*           , 

1  -  0.52  A 

t 

4 

* 
« 
* 

9 

<0 

r  •   .....  i 

* 

* 

* 

» '  ! 

I 

1     1 

:■:•::  :•*,-. 

I 1 

/N5TRO\ 

Vic  -  10  VDC 

/ 
2 

J 

* 
< 

i           1 

!'•'             1 

\.:.,::,...:..i 

0  <  V 

S  <  25  VDC 

0  <  IS      <  0.4  A 

l 

• 
9 

1 1 

INSTRON  *  2 

Figure   C-l  :     Load  Cell  Sets  connected  in  parallel. 
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Vic  -   10  VDC 
LOAD  CELLS 

/ 
2 
J 
4 

1 

1  =  0.257  A 

— i 

X 

* 
* 

9 

* 

% 

< 

........  ., 

INSTRON  *V 

/N5TRON  *  2 

Vic  =   10  VDC 

5WLTCH    y^   OFF 

Vrf 

0  < 
0  < 

V9 

Is 

< 

< 

25  Vt 

0.4  / 

)C 

Figure   C-2  :     Load  Cell  Sets  connected  in  series. 
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more  than  two  (2)  LC  sets  can  be  powered  by  one  (I)  HP-62I6B  power 
supply  as  shown  below: 

Switch  ON  :29i=,    (  1/R)  =  1/Rtot  =  9*0/350)  =  0.257  (l/fl) 

Rtot  =  38.888  0 

he =  vicREC/Rtot =  ]  °v/38  888  Q  "  °-257  A 

Switch  OFF:Vs  >  V1lc  +  V2]c  =  2  *  (VI  ]c)  =  2  *  (R  tot  *  l]c) 
V1lc  +  V2]c  =  2*38.888*0.257  V  =  19.989  V 
Hence  MAX  Vs  =  25  V  >  1 9.989  V 

However  this  conclusion  does  not  take  into  account  the  possible  failure  of 
one  LC  (short-circuit).  If  this  is  the  case  the  rest  of  the  LCs  will  be 
possibly  ruined  and  consequently  it   is  not  worthy  of  taking  the  risk. 

Therefore  it  is  decided  to  use  four  (4)  power  supplies,  one  for 
every  LC  set. 

b         LVDT  case 

According  to  the  specifications  of  the  LVDT  the  recommended 

excitation  voltage  is  24V  and  the  highest  output  resistance  106Q. 

Obviously  in  this  case  the  series  solution  does  not  work  since 
the  maximum  capacity  of  the  power  supply  is  only  25V  and  therefore  only 
parallel  connection  schemes  need  to  be  examined.  Under  the  specification 
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conditions  the  loop  is  limited  by  the  400mA  current  and  hence  by  a  total 
more  than  two  (2)  LC  sets  can  be  powered  by  one  (1)  HP-6216B  power 
resistance  of  60  Q  as  shown  below: 


MAX  ls  =  VsREC/Rtot  or  Rtot  =  VSREC/MAX  ls  =  24V/0.4A  =  600 


Hence  one  possible  yet  economic  situation  is  the  one  shown  in 
Figure  C-3  where  the  four  (4)  LVDTs  are  connected  in  parallel  to  the 
power  supply  yieling  a  total  resistance  of  600.  According  to  this  situation 
each  LVDT  amounts  for  a  resistance  of  240  0  «  106  Q  and  a  current  of 

100mA.    This  fact  means    a  voltage  differential  of  24V    which  is  well 
within  the  capabilities  of  the  HP-6216B  power  supply. 

Therefore  given  that  the  LVDTs'  reliability  is  very  high  and 
that  their  wiring  connections  will  be  spatially  that  far  appart,  that  the 
short  circuit  probability  is  very  remote,  one  (1)  power  supply  is  adequate 
for  all  the  LVDTs  that  will  be  externally  mounted  to  the  fixed  cross-heads 
of  the  INSTRON  machines. 
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Imax  =  0.4  A 

!  ■                                                         1 

i  -  l  vur  *  t        i 

j;"''    :  /•i/or*r)-         ;-| 

m 

Vs  <  25  VDC 

IS  <  0.4    A 

j.. ;.;..;  L.  YuT.    rJT 

I  ■     f'  \srtT  *",-(■             

24    VDC 

0  < 
0  < 

Figure    C-3  :     LVDTs  connected  in  parallel. 
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2.         IMPLEMENTATION  OF   THE  OVERALL  EXPERIMENTAL  SET-UP 

The  necessary  procedures  for  the  implementation  of  this  task  follow 
in  a  check  list  format.  However  in  case  clarifications  are  needed  the 
manuals  have  to  be  consulted  accordingly  [Refs.  14,  15,  16  &  17]. 


a.  Identification  of  the  appropriate  LTs  to  be  installed  in  the  moving 
cross-head  of  the  INSTRONs.  In  this  case  the  1000  lbs  LTs  had  to 
be  used  since  the  order  of  magnitude  of  the  load  at  the  end  of  the 
loading  sequence  for  the  test  of  nine  (9)  fiber  bundles 
simoultaneously  has  been  determined  to  be  of  the  450  lbs  order  of 
magnitude. 

b.  Removal  of  the  existing  100  lbs  LTs  from  the  INSTRONs" 
cross-heads,  removal  of  the  LVDTs  from  the  100  LTs  and 
installation  to  the  1000  LTs.  Installation  of  the  1000  LTs  into  the 
moving  cross-heads  and  calibration  of  the  equipment  in  accordance 
with  INSTRON  1000  manual  [Ref.  15],  (  Figure  A-4  ). 

c.  Installation  of  the  Load  Cells  (LC)/Upper  Differential  Mechanism 
System  (UDMS)  into  the  moving  cross-heads  and  of  the  Lower 
Differential  Mechanism  System  (LDMS)  into  the  lower  end  of  the 
fixed  cross-heads  as  shown  in  Figures  C-4  &  C-5. 

d.  Connection  of  the  LC  wirings  to  the  specifically  fabricated 
Distribution  Boards  (DB)  by  soldering  (  Figure  C-6a  ). 

e.  Connection  of  the  DBs  to  the  HP-3421A  DACUs  and  to  the 
HP-62I6BS  (excitation  voltage  of  the  LCs)  by  appropriately 
marked  wiring,  according  to  the  correspodence  given  in  Table 
VI. 

f.  Mounting  of  the  LVDTs  to  the  INSTRONs  fixed  cross-heads  by 
utilizing     the     specif  icaly     designed     attachement     pads     in 
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Figure   C-4  :     Load  Cells/Upper  Differential  Mechanism  system. 
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Figure    C-5  :     Lower  Differential  Mechanism  system. 
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Figure   C-6a      :     Load  Cells' Distribution  Board. 


Figure    C-6b      :     LVDTs*  Distribution  Board. 
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TABLE    VI 

LOAD  CELLS  &  LVDT  TO  HP-3421A  DACU 
WIRING  CONNECTIONS 

HP -3421 A      DACU 


SLOT     0  (LOWER  BUS) 


SLOT      1    (UPPER  BUS) 


CHANNEL  *    CONNECTED   ITEM 


CHANNEL  *    CONNECTED  ITEM 


0 

Actuator 

1 

Actuator 

2 

LC  Exctn  Vltge 

3 

Load  Cell*! 

4 

LoadCell*2 

5 

LoadCell*3 

6 

LoadCell*4 

7 

Load  Cell^S 

8 

Load  Celine 

9 

Load  Cell*? 

0 

LoadCell'8 

1 

LoadCell-r9 

2 

LVDT 

3 

LVDT  Exctn  Vltge 
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predetermined  locations  so  that  no  negative  voltage  outputs  could 
be  obtained  that  could  delay  the  data  acquisition  rate  (  Figure 
C-7) 

g.  Connection  of  the  LVDTs  wiring  to  another  specifically 
fabricated  distribution  board  (db)  by  soldering  (  Figure  C-6b  ). 

h.  Connection  of  the  db  to  the  HP-6216B  (excitation  voltage  of  the 
LVDTs)  and  to  the  HP-3421A  DACUs  by  appropriately  marked 
wiring  according  to  the  correspondence  given  in  Table  VI. 

i.  Connection  of  the  HP-85  desk  top  computer  to  the  HP-3421A 
DACUs  through  the  HP-IL  for  the  monitoring  of  the  loading 
sequence. 
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